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ABSTRACT 


Most  physical  systems  are  nonlinear  to  some  extent.  One  variety 
of  nonlinear  system  uses  a  relay  or  other  on-off  device  as  a  control 
e lement . 

Since  in  the  process  industries  or  relay  servo  systems,  one  often 
encounters  a  type  of  time  delay  called  transportation  lag  or  dead  time. 

It  is  highly  desirable  to  investigate  the  effects  of  the  transportation 
lag  on  the  process  or  relay  servo  system. 

This  paper  is  to  demonstrate  the  effects  of  the  transportation  lag 
on  the  relay  servo  system  and  show  how  to  eliminate  these  effects. 

Experimental  results  indicate  that  the  transportation  lag  causes  a 
continuous  oscillation  with  a  steady  state  over-  and  under-shoot,  the 
response  can  never  be  settled,  and  also  it  leads  to  a  longer  rise  time. 

In  addition  to  those,  the  transportation  lag  results  in  a  terminal  limit 
cycle.  Though  the  tachometer  feedback  has  been  used  to  compensate  the 
system,  this  oscillation  can  never  be  eliminated  completely  when  some  delay 
exists.  Therefore,  a  five  percent  overshoot  and  undershoot  was  considered 
acceptable . 

The  author  wishes  to  express  his  appreciation  for  the  advice,  assist¬ 
ance  and  encouragement  given  by  Dr.  Richard  Carl  Dorf,  Professor  of  Electri¬ 
cal  Engineering,  United  States  Naval  Postgraduate  School,  throughout  this 
investigation.  He  also  owes  a  debt  of  gratitude  to  Dr.  John  R.  Ward, 
Professor  of  Electrical  Engineering  Department,  U.  S.  Naval  Postgraduate 
School,  for  his  assistance  during  the  initial  stages  of  the  work. 
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1.  Introduction. 


Systems  which  utilize  feedback  for  control  purpose  have  become  es¬ 
sential  elements  in  modern  technology.  They  range  all  the  way  from 
simple  toys  to  the  most  complex  automatic  factories  and  production  equip¬ 
ment.  Feedback  control  is  the  guidance  that  makes  modern  automation  pos¬ 
sible,  and  as  such  is  in  a  large  measure  responsible  for  the  ever-increas¬ 
ing  productivity  and  rising  standard  of  living  of  man. 

It  is  noted  that  feedback  control  systems  are  seldom  truly  linear, 
but  rather  they  are  subject  to  the  introduction  of  a  considerable  number 
of  nonlinearities.  Some  nonlinearities  may  be  introduced  intentionally  by 
the  system  designer.  An  example  of  this  might  be  the  selection  of  a  re¬ 
lay  servomechanism  with  its  attendent  nonlinearity.  On  the  other  hand 
there  may  be  no  way  of  avoiding  a  nonlinearity  in  the  system  design.  An 
example  of  this  would  be  the  mandatory  inclusion  of  a  gear  train  with  its 
attendent  backlash.  A  nonlinearity  which  frequently  can  neither  be  avoided 
nor  neglected  is  tfansportat ion  lag. 

In  some  control  systems,  especially  in  the  field  of  process  control, 

a  non-eligible  propagation  time  is  encountered.  That  is,  there  is  a  definite 

lapse  of  time  between  the  initiation  of  a  change  in  the  signal  at  one  point 

in  the  system  and  the  appearance  at  some  other  point  of  any  effect  of  this 

change.  But  be  sure  that  this  kind  of  lag  is  not  to  be  confused  with  the 

^  2 

type  of  lag  introduced  by  factors  of  the  form  1/(TS  +  1)  or  1/ vT  S  +2  S  TS 
+1)  in  the  transfer  function  relating  the  two  signals.  Factors  such  as  these 
in  the  transfer  function  of  a  component  do  reveal  a  tendency  of  the  com¬ 
ponent  to  retard  the  signal,  in  a  sense,  by  introducing  phase  lag  at  all 
frequencies,  but  they  do  not  result  in  a  finite  dead  time  which  must  elapse 
before  any  effect  of  a  change  in  input  can  be  noted  in  the  output  of  the 
component . 
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Dead  time  of  appreciable  magnitude  can  be  observed,  for  example  in 
long  pipes  carrying  fluid  under  pressure,  A  change  in  pressure  at  one  end 
of  such  a  pipe  is  transmitted  with  finite  velocity  to  the  other  end,  pos¬ 
sibly  but  not  necessarily  being  distorted  in  transmission.  Another  type 
of  system  in  which  dead  time  is  often  encountered  is  shown  in  Fig,  1-1, 

In  such  a  system  two  different  fluids  are  mixed  at  the  mixing  point,  and 
the  relative  concentration  is  read  at  the  metering  point.  If  the  meter¬ 
ing  point  is  located  at  an  appreciable  distance  from  the  mixing  point, 
the  dead  time  between  the  establishing  of  a  change  in  relative  concentra¬ 
tion  and  the  detection  of  the  change  may  be  appreciable. 

Another  example  is  illustrated  below.  In  Fig,  1-2  the  shaded  area 
represents  a  slug  of  water  at  position  1,  the  outlet  of  the  steam  nozzle. 
Position  2  represents  where  the  temperature  sensing  element  will  be  placed. 
Let  the  temperature  of  the  shaded  area  water  slug  at  the  steam  nozzle  be 
quickly  changed  from  Q,  to  ,  Before  the  change  can  be  measured,  the 
slug  of  water  must  travel  from  the  steam  nozzle  to  position  2  where  the 
sensing  element  would  be  located. 

The  effect  of  the  delay  caused  by  the  slug  traveling  down  the  pipe 
can  be  shown  by  locating  sensing  elements  at  position  1  and  2  as  shown  in 
Fig.  1-3,  If  the  time  constants  of  the  sensing  element  at  position  1  is 
taken  as  the  input  and  the  element  at  position  2  as  the  output,  Fig.  1-4 
becomes  the  reaction  curve  for  the  process. 

This  type  of  reaction  curve  is  typical  of  a  dead  time  lag  (or  trans¬ 
portation  lag)  and  may  be  defined  by  the  time  lag  T.  The  magnitude  of 
the  dead  time  lag  will  depend  on  the  water  velocity  and  the  distance  be¬ 
tween  position  1  and  2,  The  lag  time  may  be  calculated  as  follows: 

T(min.)  =  d(ft)  /  V(ft/min) 
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Fig,  1-1 


Fig.  1-2 


-  input 

-  output 


t  ime 


Fig.  1-4 
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where*  T  =  dead  time  lag  (min) 
d  =  length  (ft)  # 

V  =  velocity  of  flowing  medium  (ft/min) 

Dead  time  lag  is  very  common  in  the  process  industries*  Because 
of  the  complete  lack  of  the  response  during  the  dead  time*  this  type  of 
lag  is  very  detrimental  to  control.  An  examination  of  the  equation  de¬ 
fining  the  dead  time  lag  shows  that  increasing  the  length  or  decreasing  the 
velocity  will  increase  the  lag  time*  and  we  shall  see  later*  make  the 
problem  of  control  more  difficult. 

To  determine  the  frequency  response  of  the  dead  time  lag*  a  sinusoidal 
wave  having  an  amplitude  of  IQVF  will  be  generated  around  a  bias  tempera¬ 
ture  0^.  At  some  low  frequency  the  curve  as  recorded  by  the  sensing  ele¬ 
ments  at  position  1  and  2  would  appear  as  shown  in  Fig.  l-5(a). 

Fig.  l~5(b)  and  5(c)  show  similar  results  for  two  successively  higher 
frequencies.  As  can  be  seen,  as  the  frequency  increases  the  period  decreas¬ 
es.  However*  the  lag  time  which  is  dependent  only  on  the  velocity  and 
length  remains  constant.  Thus  the  magnitude  of  the  phase  angle  will  in¬ 
crease  proportionally  to  each  decrease  in  period  (increase  in  frequency). 

The  phase  angle  may  be  calculated  as  follows: 

f>  (degree)  =  periodlScycle)  X  360  (ddgree/cycle) 
or  0  (degree)  =  -  T(min)  X  f(cpm)  X  360  (degree/cycle) 

From  this  equation  we  see  that  the  phase  angle  for  the  dead  time  lag 
is  directly  proportional  to  the  frequency.  The  phase  angle  approaches 
zero  degrees  as  the  frequency  approaches  zero  and  negative  infinity  as 
the  frequency  approaches  infinity.  This  will  be  discussed  in  jp  6. 
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2.  Relay  servo  system, 

2,1  General  Description. 

When  economy  is  of  importance,  many  feedback  control  systems  utilize 
a  contactor  for  energizing  the  actuator.  Such  a  system  is  called  a  contac¬ 
tor  or,  alternatively,  an  on-off  control  system,  since  the  contactor  is  an 
element  whose  output  can  usually  have  only  two  values.  In  some  cases 
these  correspond  to  a  fixed  positive  action  and  no  actuator  action.  This 
usually  is  the  case  with  thermostatic  controls,  for  example.  In  other 
cases  the  two  values  correspond  to  maximum  positive  and  negative  actuator 
action. 

In  Fig.  2-1  it  is  assumed  that  N  is  a  relay  or  switch  which  has  an 
output  +1  or  -I,  depending  on  the  sign  of  E  =  R  -  b.  A  positive  output  cor- 
responds  to  a  positive  sign  of  E  and  a  negative  output  corresponds  to  a 
negative  sign  of  E.  When  the  switch  is  a  perfect  switch  this  can  be  de* 
noted  symbolically  by  the  equation 

V  Sign  E  (2.1) 


Fig.  2-1  Block  diagram  cf  the  contactor  servomechanism 
Since  the  contactor  is  frequently  an  electrically  controlled  or  equivalent 
element,  there  may  be  a  delay  in  the  initiation  of  switching  action  as 
well  as  a  finite  time  interval  during  which  switching  action  takes  place. 
The  contactor  action  ordinarily  applies  full  control  force  or  torque  to 
the  system.  Since  the  output  of  N  is  therefore  not  necessarily  unity,  the 
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actual  value  is  included  in  the  gain  constant  of  the  following  element. 

The  block  G  in  Fig.  2-1  will  be  assumed  to  be  a  linear  servomotor 
with  viscous  damping  and  no  load  torque  except  that  due  to  inertia.  The 
limiting  case  of  zero  viscous  damping  is  of  interest  and  will  also  be 
considered.  Including  the  gain  constant  of  the  switch,  the  motor  transfer 
function  is 


G  = 


To 


TS(  3  f~  </€,„) 


(2.2) 


where  T  =  KE  is  the  stall  torque  of  the  motor  and  J  is  the  total  inertia 
o  o 

of  motor  and  load.  In  the  absence  of  damping,  1  t'£nj  is  set  to  zero  in  this 
equation.  It  should  be  emphasized  that  since  N  is  nonlinear  element,  there 
is  no  particular  reason  to  assume  that  the  feedback  element  is  linear.  In 
particular,  the  feedback  signal  could  be  any  algebraic  function  of  position 
and  its  derivatives.  However,  for  the  method  of  analysis  used  here  the 
feedback  function  can  depend  only  on  the  position  and  velocity  and  must 
not  on  the  derivative  of  the  velocity  (acceleration).  The  symbol  H,  there¬ 
fore,  no  longer  represents  a  linear  transfer  function,  but  a  possibility 
nonlinear  feedback  function  involving  the  output  position  and  velocity. 
Symbolically  this  will  be  written  as 


b  =  f(C,C) 


(2.3) 


where  C  =  dC/dt 
2.2  Relay  servos 

The  system  behavior  can  be  prescribed  by  combining  the  equations: 


o 

C 

+ 

/ 

* 

c 

c 

+- 

/ 

7~ 

c 

+ 


To 

J 

7Z 

J 


R  -f(  c,6)  >  o 

(2.4) 

R  -  f(cJ  c)  <c 


Between  switching  times  the  response  is  described  by  a  second  order 
differential  equation.  For  this  reason  such  a  system  is  called  a  second- 
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order  contactor  servo  or  feedback  system.  The  two  equations  are  inte- 
grable  between  switching  times  and  the  entire  solution  can  be  obtained 
by  matching  final  and  initial  conditions  of  the  two  sets  of  solutions  at 
the  proper  switching  time.  The  solution  is  therefore  easy  to  obtain  in 
principle .  The  details,  however,  are  so  complex  that  another  method  of 
solution  is  desirable. 

The  analysis  of  the  system  can  be  simplified  by  using  the  fact  that 
it  is  possible  to  plot  the  response  of  the  system  as  a  trajectory  in  the 
phase  plane,  the  trajectory  being  composed  of  a  series  of  arcs  or  curves 
which  have  a  characteristic  shape  dependent  upon  the  properties  of  the 
system.  The  behavior  of  a  physical  system  where  only  one  coordinate  is  of 
interest  is  usually  specified  by  determining  this  coordinate  as  a  function 
of  time,  C  =  fj  (t).  This  function  depends  upon  the  differential  equation 
describing  the  system  and  upon  a  set  of  initial  conditions.  For  some 
purpose,  an  equal  suitable  way  of  describing  the  behavior  of  the  system 
is  by  means  of  a  functional  relationship  between  position  and  velocity: 


(2.5) 


When  such  a  relation  is  known,  the  time  dependence  of  C,  if  required, 
can  be  found  directly  by  integrating  the  equation 


t/g 

cit 


c 


7 C  (c.) 


Or 


rC 


</_JL 

fCv 


(2.6) 


Equation  2.5  can  be  represented  graphically  with  position  and  velocity 
taken  as  rectangular  cartesian  coordinates;  position  is  usually  plotted 
as  abscissa  and  velocity  as  ordinate.  The  behavior  of  the  system  is  des¬ 
cribed  by  a  plot  of  all  successive  pairs  of  values  of  position  and  velo¬ 
city.  The  plane  on  which  this  plot  is  made  is  called  the  phase  plane,  and 
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the  graph  of  Equation  2.5  is  called  the  phase-plane  trajectory. 

For  the  system  to  be  considered  in  this  paper,  the  differential 
equations  of  interest  will  not  contain  time  explicitly.  For  the 
second-order  system  this  implies  that  the  acceleration  can  be  written  as 

c  =F(c,  c) 

The  importance  of  this  concept  here  results  from  the  fact  that  for  such 

systems,  the  behavior  and  therefore  the  trajectory  are  uniquely  specified  by 

# 

the  initial  values  of  C  and  C.  This  means  that  there  is  an  only  tra¬ 
jectory  through  a  given  (nonsingular)  point  in  the  phase  plane.  In  this 
paper  it  is  restricted  to  the  systems  described  by  second  order  differential 
equations  and  hence  by  phase  plane  trajectory  in  the  two  dimensional  C-C 
phase-plane. 

The  slope  of  the  trajectory  in  the  phase  plane  can  be  written  as 

dc 

dc  d  t 

dc  dc  (2,7) 

dt 

According  to  this  equation  the  slope  is  equal  to  the  quotient  of  ac- 

«» 

celeration  and  velocity.  In  the  upper  half  plane,  C  is  positive  and 
nonzero,  and  therefore  the  slope  is  nonifinite  and  either  positive  nega¬ 
tive,  depending  upon  whether  the  acceleration  is  positive  or  negative. 
Similarly,  in  the  lower  half  plane  the  slope  is  non-infinite  and  either 
positive  or  negative,  depending  on  whether  the  acceleration  is  negative 

or  positive.  For  a  point  on  the  horizontal  axis  the  velocity  is  zero, 

* 

and  therefore  unless  the  accelerations  also  zero  the  slope  dC/dC  is  in¬ 
finite.  In  general,  then,  trajectories  cross  the  horizontal  axis  verti¬ 
cally.  Direction  of  motion  along  the  trajectory  can  always  be  determined 

e 

by  the  fact  that  for  C  positive  C  must  become  more  positive  with  time  and 
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for  C  negative  C  must  become  more  negative  with  time.  Motion  along 
the  trajectory  is  therefore  to  the  right  in  the  upper-half  plane  and  to 
the  left  in  the  lower-half  plane.  The  horizontal  component  of  the 
trajectory  velocity  is  equal  to  the  system  velocity  C  =  dC/dt  and 
the  vertical  component  of  the  trajectory  velocity  is  equal  to  the  system 

c 

acceleration  C  =  dC/dt.  These  considerations  make  it  clear  that  a  trajec 

tory  can  not  end  except  on  the  horizontal  axis.  If  there  are  points  on 

*  •  •> 

the  horizontal  axis  where  both  C  and  C  are  zero  they  are  called 
singular  points,  and  the  value  of  C  at  such  a  point  corresponds  to  an 
equilibrium  condition  of  the  system.  Trajectories  can  emerge  from  or 
return  to  some  of  these  singular  points  with  a  non- infinite  slope,  since 
the  slope  dC/dC  is  undefined  at  this  point  except  through  other  considera 
tions, 

2.3  Phase  plane  trajectory. 

Some  of  the  ideas  in  ^2,2  will  now  be  illustrated  with  the  system 
shown  in  Fig.  2.1  and  defined  by  Equation  2.2  through  2.4,  The  Equation 
2,4,  which  are  second-order  differential  equations,  can  be  reduced  to 
pairs  if  first-order  equations  by  substitution  of  Equation  2.6 

k=  C  =  r/C 

which  gives 


Cft 


i/  f  /  i  /  /  7c 


(2.8) 


an  alternative  expression  for  velocity,  dV/dt  =  V(dV/dC),  can  be  used 
with  Equation  2,8  to  give  the  following  equation: 


1/ 


d  l 

dc 


V 


% 


-  J 

This  equation  can  now  be  integrated  with  the  result 

c  7V 


(2.9) 


c-c.=  v  LfO  +£%)] 
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Evidently  the  constant  of  the  integration  is  the  value  of  C  at 


which  V  =  0.  This  equation  can  be  put  in  dimensionless  form  by  the 
substitutions 

J 


1/  =  V 

C„To 


C 


?jr„ 


c 


which  gives 

z  -z0 


-[  V  *  ^  d  +  I 7)J 


(2.11) 


(2.12) 


Typical  trajectories  for  particular  values  of  art  shown  in  Fig,  2,2. 
Through  use  of  the  dimensionless  variables  v  and  c$  all  systems  of  this 
type  have  identical  trajectories. 

In  the  absence  of  damping,  the  trajectories  are  parabolas.  If 
is  set  to  zero,  integration  of  Equation  2.8  gives  the  result 


c  -  Co  =± 


X 


[/* 


2 


(2.13) 


In  this  case  the  substitution  2.11  can  not  be  used.  For  the  case  of 
finite  damping  it  is  sometimes  convenient  to  have  an  expression  for  the 
trajectory  which  shows  the  deviation  from  the  parabolic  shape  more 
directly  than  Equation  2,10  or  2.12.  This  can  be  obtained  by  use  of  the 


power  series  expansion  for  log  (1  +  x): 

To 


c_c  ==tJLJZ 

W  V-o  *4  -7-  £ 


[It  t 


TmTo 


_£  /«7V 

4  *  Z,„  To 


A -7 


or 


c  -  c0  =  t 


o* 


<4 


o  2 

~  (PJt 


(2.14) 


(2.15) 


These  equations  are  particularly  useful  when  an  approximate  result  is 
desired.  For  accurate  computations  Equation  2.10  or  2.12  is  more  use¬ 
ful,  since  the  series  given  above  converges  rather  slowly  except  for 
small  values  of  v,  and  indeed  diverges  for  jvj  ^  1. 

Reference  to  Fig.  2.1  and  Equations  2.1  and  2,3  shows  that  the 
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Fig.  2.2,  Phase  plane  trajectory  for  the  second-order 
contactor  servomechanism. 

condition  governing  switching  from  one  trajectory  family  to  the  other  is 
given  by  the  equation 

/?  -  f  (c,  C)  =  o  (2.16) 

For  the  present  it  will  be  convenient  to  restrict  attention  to  step  in¬ 
puts  with  initial  velocity  zero.  In  this  case  a  simple  translation  of 
the  origin  can  be  used  so  that  the  initial  conditions  are  c  =  -c  ,  V  =  0, 
and  the  final  desired  steady-state  conditions  are  R  =  C  =  V  »  0.  Setting 
R  =  0  in  Equation  2.16  gives  the  equation 

(C*  £  )  =  O  (2.17) 

The  graphical  plot  of  this  equation  separates  the  phase  plane  into  two 
parts,  in  one  of  which  the  system  follows  trajectories  corresponding  to 
V  -  +1,  and  in  the  other  of  which  the  system  follow  trajectories  corres¬ 
ponding  to  V  =  -1,  this  is  illustrated  in  Fig.  2.3. 

For  certain  functions  f(c,c)  an  anomalous  situation  arises.  Consider, 
for  example,  the  point  P^  in  Fig.  2,3,  which  will  supposed  to  be  reached 
by  motion  along  the  trajectory  marked  (1).  When  the  switching  line  is 
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V 


Fig.  2-3.  Switching  line  for  contactor  servomechanism. 
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reached*  the  system  will  normally  transfer  to  the  trajectory  of  the 
other  family  which  pass  through  P  ,  this  is  the  trajectory  (2)  in  Fig. 

2.3.  However*  when  the  point  is  reached*  for  example  along  trajectory 
(3),  a  paradoxal  situation  arises.  Although  switching  should  take  place 
at  this  time,  there  is  no  trajectory  of  the  other  family  along  which  the 
system  can  move,  since  P^  is  an  end  point  rather  than  a  starting  point  for 
trajectory  (4).  In  terms  of  real  physical  systems  the  paradox  can  be 
resolved  by  recognizing  that  a  small  delaying  switching  is  unavoidable. 
This  permits  the  system  to  overshoot  the  switching  line  so  that  a  small 
section  of  allowable  trajectory  is  available  in  other  family.  When  this 
trajectory  reaches  the  switching  line  another  slight  delay  allows  over¬ 
shoot  back  into  the  left-half  plane.  The  resultant  trajectory  consists 
of  a  series  of  short  arcs  in  which  the  system  oscillates  about  the  switch¬ 
ing  line  as  it  approaches  the  equilibrium  point  at  the  origin.  In  terms 
of  the  ideal  system,  one  can  say  that  the  trajectory  follows  the  switching 
line  to  the  origin  accomplished  by  an  infinite  number  of  switching  which 
provide,  on  the  average,  the  proper  force  to  give  an  acceleration  consist¬ 
ent  with  the  location  in  the  phase  plane  and  the  slope  of  the  switching 
line . 

Using  the  above  considerations  some  interesting  general  conclusions 
can  be  drawn: 

(1)  For  no  viscous  damping  (  —  oa)  and  no  rate  feedback 

[•/Yc  c) ^c]  the  trajectories  are  parabolas  separated  by  the  vertical  line 

c  =  0.  The  system  is  neutrally  stable  in  the  sense  that  for  any  amplitude 
of  step  input  the  oscillation  will  persist  indefinitely,  as  shown  in  Fig, 
2.4a. 

(2)  For  viscous  damping  and  no  rate  feedback  j  f  the 
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c 


(a) 

c 


J 

c 


Fi&°  2'4*  Typical  phase  plane  trajectories 
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trajectory  are  the  "asymmetric  parabolas"  of  Fig.  2.2,  separated  by  the 
verticle  line  c  =  0.  The  system  exhibits  a  damping  or  decrement  of  the 
oscillations  produced  by  a  step  input.  The  per-unit  decrement  is  a  func¬ 
tion  of  amplitude  and  approaches  zero  as  the  amplitude  approaches  zero. 

In  other  words,  for  large  input  the  initial  per-unit  decrement  is  high,  but 
decreases  to  zero  as  the  amplitude  of  the  oscillation  approaches  zero. 

This  can  be  observed  qualitatively  in  Fig.  2.4b. 

(3)  With  linear  rate  feedback,  the  feedback  function  is  ffc* 

If  there  is  no  viscous  damping  (  ,  then  the  trajectories  are 

parabolas  separated  by  the  straight  line 

C  =  -  -fi-  C  (2.18) 

Again  the  system  shows  a  per-unit  decrement  or  damping  which  varies  with 
amplitude.  However,  in  this  case  the  damping  is  small  for  large  amplitude 
and  as  the  amplitude  decreases  the  damping  increases  until  the  final  tra¬ 
jectory  into  the  origin  lies  along  the  line 


with  no  further  overshoots.  An  illustration  of  this  is  shown  in  2.4c. 

The  above  conclusions  shows  that  no  simple  linear  feedback  function 
will  give  optimum  response,  regardless  of  the  system's  viscous  damping. 

If  rate  feedback  is  adjusted  to  give  no  overshoot  for  the  largest  expected 
input  then  the  system  moves  to  the  origin  slowly  for  small  inputs.  If  the 
rate  feedback  is  adjusted  to  give  faster  response  for  small  input,  then 
overshoot  and  oscillations  will  follow  application  of  a  large  step  input. 

It  is  natural  to  ask  whether  any  feedback  function  can  be  found  to 
give  optimum  response  for  all  inputs.  This  turns  out  to  be  a  difficult 
question.  Since  the  system  is  nonlinear  and  does  not  have  an  easily  deter¬ 
mined  response  for  arbitary  inputs.  However,  a  definite  answer  can  be  found 
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for  the  case  of  step  inputs  by  following  chain  of  reasoning.  For  most 
rapid  response  without  overshoot  it  is  necessary  to  apply  full  positive 
acceleration  torque  until  that  instant  at  which  full  negative  torque  will 
bring  the  system  to  rest  at  the  origin  in  the  phase  plane.  This  is  il¬ 
lustrated  in  Fig.  2,5  for  a  particular  initial  starting  point.  The  tra¬ 
jectory  shown  passing  through  the  origin  can  be  designated  as  the  homing 
trajectory.  In  order  to  achieve  switching  on  this  trajectory  for  all  in¬ 
puts,  it  is  necessary  that  feedback  function  set  equal  to  zero, 

j  (<z,  a)  —  O  (2.19) 

being  the  equation  of  this  homing  trajectory.  In  this  case,  when  any 
trajectory  crosses  this  line  switching  will  occur,  so  that  the  applica¬ 
tion  of  full  reverse  torque  will  just  bring  the  system  to  reset  at  the 
origin. 


Fig.  2.5  Switching  point  for  optimum  step  response 
Since  the  trajectory  of  Equation  2.8  is  roughly  parabolic,  a  suit¬ 
able  form  of  the  feedback  function  is 

'f  Cc'  v)  ^  k  ( C  t  \  )=  C  ~+  V  ji/| 
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Strictly  speaking,  this  equation  should  be  modified  to  match  the  appropri- 
ate  "asmmetric  parabola”,  although  for  practical  purposes  the  refinement  is 
not  necessary.  In  other  words,  if  the  position  feedback  is  linear  and  the 
rate  feedback  roughly  quadratic  (except  for  sign)  then  it  is  possible  to 
obtain  optimum  response  for  this  particular  type  of  input.  It  will  be 
seen  that  switching  delays  give  rise  to  oscillatory  motion  of  small  ampli¬ 
tude  unless  a  linear  velocity-feedback  term  is  introduced.  For  this  reason 
it  is  usually  desirable  to  introduce  a  small  amount  of  such  feedback  (in 
addition  to  the  quadratic  feedback)  in  order  to  eliminate  or  minimize  this 
oscillation.  This  does  not  significantly  affect  the  results  for  large 
amplitudes . 

It  is  quite  important  to  remember  that  this  optimization  does  not 
necessarily  give  optimum  response  for  arbitary  inputs.  Further  considera¬ 
tion  is  beyond  the  intended  scope  of  this  paper. 

2.4  On-off  control. 

There  probably  are  more  on-off  discontinuous  controllers  in  the  world 
than  the  number  of  all  other  servomechanisms  and  regulators  put  together. 
Hot-water  heaters,  ovens,  refrigerators,  furnaces,  electric  irons,  battery 
chargers,  liquid-level  controls,  and  compressed-air  pressure  regulators 
are  all  commonly  of  this  type. 

A  relay  system  is  one  which  contains  a  decision  function.  It  is 
commonly  in  the  error  amplifier  and  in  its  most  elementary  form  is  an  on- 
off  switch.  If  the  positive  error  exceeds  a  certain  amount,  the  switch  is 
turned  from  on  to  off,  and  if  the  negative  error  exceeds  a  certain  amount, 
the  swtich  is  turned  from  off  to  on.  When  the  switch  is  on,  a  predetermined 
constant  corrective  force  is  applied  into  the  rest  of  the  system  to  bring 
the  output  back  into  correspondence  with  the  desired  value.  When  the  switch 
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is  off,  this  force  is  removed.  The  differential  equation  for  each  is  the 
open-loop  equation,  with  one  term,  the  corrective  force,  having  two  differ¬ 
ent  values  for  the  two  periods.  Each  switch  operation  starts  a  new  period 
for  which  the  initial  values  are  the  final  values  of  the  previous  period. 

One  characteristic  of  a  relay  is  that  it  delivers  an  output  having  one 
of  only  two  values,  independent  of  the  continuous  changes  in  the  input, 
except  they  specify  the  switching  instants.  It  also  has  hysteresis,  the 
phenomenon  that  the  input  values  for  switching  on  and  for  switching  off 
are  not  exactly  equal.  The  actual  time  of  motion  of  the  armature  will  be 
considered  negligible,  but  the  system  may  contain  the  dead  time  of  flow  of 
material  through  a  long  pipe,  equivalent  to  the  time  of  transmission  of  an 
electric  signal  on  a  long  distoritonless  line. 

The  term  relay  is  used  to  apply  to  any  device  having  at  least  the  first 
of  these  characteristics.  In  this  paper,  we  are  interested  primarily  in 
discontinuous  control  of  continuous  processes,  so  that  the  relay  is  only 
one  part  of  the  system. 

The  on-off  control,  or  two  position  control,  is  undoubtedly  the  most 
widely  used  type  of  control  for  both  industrial  and  domestic  service.  It 
is  the  kind  of  control  generally  employed  on  home  heating  systems  and 
domestic  water  heaters. 

Two  position  control  is  a  position  type  of  controller  action  in  which 
the  manupulated  variable  is  quickly  changed  to  either  a  maximum  minimum 
value  depending  upon  whether  the  controlled  variable  is  greater  or  less 
than  the  set  point.  The  minimum  value  of  the  manupulated  variable  is  usually 
zero  (off).  This  mode  of  control  is  illustrated  by  the  electric  level- 
control  in  Fig.  2,6.  A  float  in  the  vessel  operates  an  electric  switch 
which  controls  power  to  a  solenoid  valve.  When  the  liquid  level  rises, 
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Fig.  2-6.  Two  position  control 


m 


Fig.  2-7 

rwo-position  control  without 
differential  gap 


m 


Fig.  2-8 

Two-position  control  with 
differential  gap 
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the  switch  contacts  are  closed,  the  solenoid  valve  closes,  and  the  flow  is 
cut  off.  When  the  liquid  level  falls,  the  switch  contacts  are  opened,  the 
solenoid  valve  opens,  and  the  inflow  resumes.  If  the  float  lever  has  no 
bearing  friction  and  the  electrical  contact  draw  no  arc,  the  action  is  sharp 
or  "knife1 edge”  as  shown  in  Fig.  2-7. 

The  equation  for  the  two-position  control  are 

m  =  when  e  P'"  0 

m  =  M  when  e  0 
o 

where  m  =  manupulated  variable 

=  maximum  value  of  manupulated  variable  (on) 

=  minimum  value  of  manupulated  variable  (off) 
e  =  deviation 

Thus  two-position  control  must  be  described  by  two  equations,  each  apply¬ 
ing  in  a  certain  region  of  deviation. 

A  differential  gap  in  two-position  control  causes  the  manupulated  vari¬ 
able  to  maintain  its  previous  value  until  the  controlled  variable  has  moved 
slightly  beyond  the  set  point.  In  actual  cperation  it  is  the  same  as  hystere¬ 
sis,  as  may  be  seen  in  Fig.  2-8. 

A  differential  gap  is  caused  in  the  two-position  controller  of  Fig.  2.6. 
If  small  static  friction  exists  at  the  bearing  on  the  float  arm,  the  liquid 
level  must  rise  slightly  above  the  desire  value  to  create  sufficient  buoyant 
force  to  overcome  friction  when  the  level  is  rising.  Also,  the  liquid  level 
must  fall  slightly  below  the  desired  value  when  the  level  is  falling  so  that 
the  weight  force  may  overcome  the  friction.  This  kind  of  differential  gap 
may  be  caused  by  unintentional  friction  and  lost  motion. 

A  differential  gap  may  be  intentional,  as  when  a  magnet  is  installed  on 
the  float  arm  in  Fig.  2.6,  causing  a  hysteresis  in  float-arm  action.  Similar 
arrangements  are  common  in  domestic  thermostats  and  are  employed  for  the 
purpose  of  preventing  rapid  operation  of  switches  and  solenoid  valves  for 
reducing  of  electrical  contacts. 
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3 .  Description  of  the  system. 


The  nonlinear  feedback  control  system  to  be  studied  is  a  relay  servo 
system  as  shown  in  Fig*  3-1. 

Block  one  contains  the  transportation  lag.  The  nature  of  this  block 
and  its  analytic  representation  will  be  discussed  in  section  6. 

Block  two  is  the  ideal  relay,  no  dead  zone  and  no  hysteresis* 

Block  three  is  the  motor  load  unit  or  plant.  It  is  assumed  to  be  a 
second  order  system*  In  order  to  have  the  transfer  function  always  low 
pass  in  nature,  we  choose  the  value  of  G(s)  as 


/ 


(3.1) 


s  (Si-t) 


Block  four  is  a  tachometer.  In  this  system  velocity  feedback  compen¬ 
sation  is  used. 

This  system,  in  addition  to  the  velocity  feedback,  will  also  use  unity 
position  feedback. 
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Fig.  3-1  Generalized  block  diagram 
of  basic  system 


4. 


Nonlinear  analysis. 


4.1  General 

In  many  cases  the  system  designer  is  forced  to  work  with  components 
which  are  not  entirely  linear.  In  order  to  linearize  the  response  of  a 
nonlinear  actuator  we  use  the  feedback.  The  basic  difficulty  with  the  non¬ 
linear  equation  is  that  it  is  not  possible  to  exhibit  a  general  solution. 
Associated  with  this  is  the  fact  that  there  is  no  complementary  function 
or  transient  solution  made  up  of  terms  with  characteristic  shapes,  there 
are,  therefore,  no  system  roots  such  as  were  defined  for  linear  systems 
with  constant  coefficients.  It  is  not  the  scope  of  this  paper  to  include 
a  comprehensive  treatment  of  nonlinear  systems. 

Many  systems  are  approximately  linear  when  the  amplitude  of  the 
system  variables  are  not  too  large  or  too  small.  For  large  enough  ampli¬ 
tude  all  physical  systems  must  become  nonlinear.  A  system  which  for  small 
amplitude  is  a  linear  unstable  system  must  become  nonlinear  when  the  oscil¬ 
latory  amplitude  exceeds  some  value.  Such  a  system  will  sometimes  reach  a 
stable  oscillatory  motion  which  is  termed  a  limit  cycle.  Although  such 
motion  is  usually  periodic,  it  need  not  be  sinusoidal. 

There  are  other  large-amplitude  nonlinearities  which  are  in  a  sense 
destablizing.  In  these  systems,  small  amplitude  response  may  be  stable 
and  satisfactory,  but  large  disturbances  may  send  the  system  into  uncontrol¬ 
led  large-amplitude  oscillations. 

There  are  also  nonlinearities  which  become  important  for  small  amplitude 
motions.  Mechanical  friction  and  backlash  are  good  examples  of  this  type. 

In  the  case  of  backlash  it  is  not  unusual  to  find  a  limit  cycle  of  small 
amplitude  which  leaves  the  system  capable  of  following  commands  or  inputs 
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with  good  accuracy  and  in  an  essentially  linear  mannes  except  for  the  small 
limit  cycle  motion  which  is  superimposed  on  the  desired  response. 

4.2  The  describing  function. 

Many  functions  involved  in  automatic  control  are  not  desirable  by 
elementary  techniques,  and  they  must  be  described  by  graphical  or  other 
methods. 

The  output  of  a  nonlinear  element  will  not  necessarily  be  sinusoidal 
even  though  the  input  is  a  pure  sinusoid,  but  it  will  generally  contain 
fundamental  and  harmonic  waves.  The  frequency  of  the  fundamental  of  the 
output  will  often  be  the  same  as  the  frequency  of  the  input.  In  this  case 
a  describing  function  can  be  used  and  the  element  may  be  treated  as  a  lin¬ 
ear  element. 

Since  most  elements  of  a  feedback  system  are  of  a  low-pass  nature,  the 
higher  harmonics  may  be  greatly  attenuated  in  the  rest  of  the  loop,  so  that 
the  input  to  the  nonlinear  element  may  be  quite  accurately  sinusoidal.  If 
this  is  the  case,  it  is  necessary  to  compute  the  transmission  around  the 
loop  for  the  harmonics.  In  particular,  for  the  nonlinear  element  we  need 
only  a  function  to  describe  the  fundamental  component  in  the  output  as  a 
function  of  the  amplitude  and  frequency  of  the  sinusoidal  input.  This  furc- 
tion  has  been  named  a  describing  function.  Since  the  describing  functions 
are  approximations  to  the  output-input  ratio  of  an  element  having  slight 
nonlinearily,  it  is  anolagous  to  the  transfer  function  of  a  linear  system, 
and  for  certain  purposes  will  be  used  in  a  similar  manner.  There  is,  of  course, 
one  important  difference.  Since  a  transfer  function  describes  a  linear  sys¬ 
tem,  it  may  be  frequency-dependent  but  it  is  not  amplitude-dependent.  On 
the  other  hand,  the  describing  function  is  usually  amplitude-dependent  and 
may  also  be  frequency-dependent.  The  describing  function  is  expressed  in 
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terms  of  the  amplitude  and  phase  of  the  fundamental  of  the  output  signal 
related  to  the  amplitude  and  phase  of  the  input  signal.  The  describing 
function  can  be  obtained  analytically  or  experimentally. 

Describing  functions  are  generally  useful  in  automatic  control  anal¬ 
ysis  if  the  neglect  of  the  harmonics  in  the  output  signal  does  not  cause 
difficulties.  They  may  be  used  with  adequate  results  when  the  element  in 
question  is  followed  by  one  or  more  elements  having  appreciable  attenua¬ 
tion  of  the  harmonics  in  question. 

4.3  Computation  of  the  describing  function. 

Consider  a  nonlinear  element  whose  input  is  sinusoidal, 


x  =  X  sin  Oft 


(4.1) 


and  whose  steady  state  output  can  be  written  as  a  Fourier  series: 


b^  sino/t  +  a^coslc^Jt  +  b^  sin  2aJt 


+ 


=  a  I2  —  ^  (a  cosk^ut  +  b  sink^t  )  +  ....(4.2) 

O ,  t-  K  K 


Evidently 


(4.3) 


b 


(4.4) 


The  describing  function  is  defined  as  the  complex  number 


N  =  Nj  +  jN2  =  Nt€ 


jf 


(4.5) 


where 


N2  =  ax/x 


or,  alternatively, 


(4.6) 
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This  definition  immediately  puts  certain  restrictions  on  the  permissible 
nonlinearities.  The  output  function  must  be  a  periodic  function  with 
period  equal  to  that  of  the  input.  According  to  this  definition,  devices 
which  produce  subharmonic  response  do  not  have  a  describing  function. 

The  form  of  Equation  4.1  implies  that  the  input  is  sinusoidal  and 
does  not  contain  d-c  bias.  If  a  nonlinearity  in  a  feedback  system  is  not 
symmetric,  then  in  general  there  will  be  such  a  bias. 

Because  we  are  only  interested  in  the  nonlinear  device  of  ideal  relay 
only  ideal  relay  will  be  discussed.  The  characteristic  of  the  ideal  relay 
is  as  shown  in  Fig.  4-1. 


+  V 


•+  Y 


+ 

-V 

°  ei 

Fig.  4-1.  Characteristic  of  ideal  relay. 

From  Fig.  4-1 


e  =  +  Y  o  art  jt 

o 

e  =  -  Y  rt  <^t<  2 re 

o 

The  describing  function  of  the  ideal  relay  can  be  obtained  very  simply 
by  carrying  the  integration  indicated  in  Equations  4.3  and  4.4  and  4.6.  And 
in  calculating  the  describing  function,  we  are  only  interested  in  the  funda¬ 
mental  terms.  Thus  we  have  the  describing  function  of  the  ideal  relay  as 


QJ/oj)= 


4  y 

7TX 


(A. 7) 
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From  the  above  results,  we  know  that  the  describing  function  of  the 
ideal  relay  is  amplitude-dependent  and  independent  of  frequency.  This 
function  is  plotted  against  X  in  Fig.  4.2. 


Fig.  4.2.  Describing  function  of  ideal  relay. 

4.4  Analysis  of  the  limit  cycle. 

4.4.1  Frequency  response  method 

Limit  cycle  angular  frequency  ,  as  predicted  by  the  frequency  re¬ 
sponse  method,  is  such  that  the  loop  phase  shift  is  360  degrees  and  this 
occurs  when  the  loop  frequency  response  function  has  no  imaginary  part. 

The  flow  diagram  for  sinusoidal  quantities  is  shown  in  Fig.  4.3;  the 
relay  is  represented  by  a  describing  function  in  cascade  with  a  time  de¬ 
lay  transfer  function  and  the  relay  output  takes  on  the  values  +  or  -  unity. 
The  condition  for 


Fig.  4.3. 

a  limit  cycle  is  therefore  met  when 


Q 


-Jos*,  td 


is  real 


</*£  (ju!±  -+  /J 

by  rationalizing  this  expression  the  condition  becomes: 


i-J & L  J 
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is  real 


Expanding  the  exponential  term  gives  the  limit  cycle  condition  as 


(Cos  -  jSin  C<JZ  td)(6jJ  +  /“U  is  real 


Equating  the  imaginary  terms  to  zero  yields  the  condition 


OJl  fee*  COi  =  / 


P 


2  7T 


(4.8) 


2n 


For  a  given  value  of  t  ,  (delay  time)  the  dimensionless  period  P  =  —— — 

d 

is  multivalued.  These  values  correspond  to  the  possible  modes  of  oscilla¬ 
tion  and  only  the  longest  period  is  stable. 

According  to  the  frequency  response  method  the  servo  loop  gain  at  limit 
cycle  frequency  is  unity  since  the  oscillation  is  exactly  regenerated  in  the 
loop.  To  determine  the  limit  cycle  amplitude  at  the  input  to  the  on-off 
element  it  is  therefore  necessary  to  solve  the  equation  obtained  by  equating 
the  loop  frequency  response  function  to  unity  at  the  angular  frequency,  that 


IS 


_  2  «  j 

n  Xu  l  JoJL  -  OJj?  J  1 


(4.9) 


where  is  the  amplitude  of  the  fundamental  component  of  the  error  signal 
when  limit  cycling.  Thus  at  limit  cycle  frequency,  the  imaginary  part  of 
the  loop  frequency  response  function  is  zero  and  by  rationalizing  the  left- 
hand  side  of  equation  4.9  and  equating  the  imaginary  part  to  zero  we  obtain 
4  j  C-Q$  td  -t  oJj. 

l 


4 

7TX, 


or 


/  ■+ 


X  =  Y  J=r  -fL  p  J  &&  Cos  tci  7 

A/.  '  J  r+  co-j  /  <4’10) 

where  x^  is  the  dimensionless  measure  of  misalignment  limit  cycle  amplitude 

and  P  is  the  time  limit  cycle  period. 

The  values  of  P  and  x  from  equations  4.8  and  4.10  are  plotted  as  shown 

L 
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in  Fig.  4, 4.  From  Fig.  4.4  we  can  predict  limit  cycle  amplitude  and  period 
in  an  on-off  system. 

Equation  4.9  can  be  rewritten  as 


G  ^  (4-11) 

Equation  4.11  may  be  analyzed  for  stability  by  Nyquies  method  or  by  root 
locus  method  if  the  describing  function  is  amplitude-dependent  only. 

When  is  amplitude-  but  not  frequency-dependent,  it  is  possible  to 
reduce  a  very  simple  criterion  for  predicting  stability  or  instability. 

The  basis  for  the  criterion  is  that  the  leftside  of  equation  4.11  is  only 
frequency-dependent  and  the  right  side  is  only  amplitude-dependent.  Margin¬ 
al  stability  as  predicted  from  the  Nyquist  criterion  and  the  describing  fun¬ 
ction  concept  requires  a  solution  of  equation  4,11,  or  in  other  words,  an 
intersection  of  these  two  curves.  At  the  point  of  intersection,  the  frequency 
parameter  on  the  one  curve  and  the  amplitude  parameter  on  the  other  give  the 
frequency  and  amplitude  of  the  corresponding  oscillation.  If  there  is  no 
intersection,  then  the  system  is  presumably  stable  or  unstable  for  all  condi¬ 
tions  of  operation.  Since  it  is  difficult  to  find  a  rigorous  definition  of 
stability  for  nonlinear  systems,  the  above  statement  must  be  regarded  as  pro¬ 
visional,  pending  more  detailed  investigation  or  more  careful  definition  of 
stability  for  urusual  systems. 

From  equation  4.11,  we  have  ■  '**  — - 

j  vj  ( -r  t) 


_ / _ 7T  X  .  _  __  . 

Qo  4-Y  /'  5  70  $  X 

if  we  put  t^  =  .2  seconds  and  by  assigning  some  values  of  zxT  ,  we  can  have 

the  following  data;  and  from  these  data  the  Nyquist  plot  can  be  plotted  as 

shown  in  Fig.  4.5. 
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Fig,  4-4  Curves  for  predicting  limit  cycle  amplitude  and 

period  in  r. n  on-off  system 
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From  the  Nyquist  plot,  a  limit  cycle  is  observed;  it  is  at  a  point  where 
,52  and  {  1*71.  (The  intersection  of  these  two  curves).  The 

intersection  may  indicate  a  stability  boundary  in  the  sense  that  excitation 
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at  level  below  that  corresponding  to  the  intersection  will  be  unstable  for 
the  system,  and  should  reach  a  limit  cycle  in  which  the  input  to  the  non¬ 
linear  element  is  roughly  sinusoidal  with  angular  frequency  OJ  ~ * 52  and 
X  =  1.71.  The  system  should  be  unstable  for  low  level  of  excitation  and 
stable  for  level  above  that  corresponding  to  the  intersection. 

4.4.2  The  root  locus  method 

Because  the  describing  function  of  an  ideal  relay  is  independent 
of  frequency,  so  that  the  root  locus  method  can  be  readily  applied.  The 
nonlinearity  introduces  no  phase  shift  but  only  a  gain  dependent  on  ampli¬ 
tude.  It  will  be  assumed  here  that  the  behavior  of  the  nonlinear  system 
can  be  approximated  through  use  of  a  root-locus  diagram  for  the  system  with 
the  nonlinearity  eliminated.  By  block  diagram  manupulation  as  follows; 


( a ) 


Fig.  4.6. 

and  algebraic  manupulation  provides  the  stability  relationship,  the 
characteristic  equation 

Gj  G2  Gd(H  +1)  =  -  1  (4-12) 
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G(S)  =  -  1 


(4.13) 


g(s)  -  2  e"st 

S(S  +  1  ) 


K'  = 


4Y 

nX 


For  t  =  .5  second, 
d 

put  Kt 


Kfc  S  =  H 


=1.5  then  we  have 


_  K  (S  +  0.67)e 
G(S)  ~  S  (S  +  1) 


- .  5S 


Equation  4.13  is  complex  and  may  be  split  into  two  equations  by  equating 
the  magnitude  and  phase  angle  on  both  sides  respectively.  Thus  we  have  the 
magnitude  equation 

| G(S)|  =  1  (4.14) 

and  the  phase  angle  equation 

/G(s)  =  ±  180°  (4.15) 

where 


where 


/G(s)  =  /g'(s)  ~  Ctrl 

/g'(s)  -Ml  =  +  180°  (4.16) 

g ' (s,  =  ■* 

V  '  S  (s  +  1) 

a  method  is  introduced  to  obtain  the  required  locus  of  equation  (4.16). 
First  it  is  to  construct  a  family  of  root  loci  (or  phase-angle  loci)  of 
the  equation 


SS  '  (s)  =  0X 

where  0^  is  a  constant  phase  angle.  Thus  a  family  of  loci  can  be  plotted 
by  various  values  of  0^,  Fig.  4.7  illustrates  a  family  of  root  loci  for  the 
transfer  function  G  t(c;).  Since  the  locus  with  phase  angle  0,  is  symmetrical 
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Fig.  4-7  Construction  of  the  Required  locus  with 

Transportation  lag 


to  that  with  phase  angle  (-0)  with  respect  to  the  real  axis,  only  the  loci 
on  the  upper  half  plane  are  shown. 

Next,  it  is  to  construct  another  family  of  phase-angle  loci  (the  name 
dead  time  lag  loci  is  suggested)  on  the  s-plane  from  the  relation 

-CO  T  =  02 

the  dead  time  lag  loci  family  is  also  plotted  in  the  same  sheet  as  shown  in 
Fig.  4.7.  Then  by  superimposing  these  two  families  of  loci  on  the  same 
s-plane,  the  point  of  intersection  of  the  corresponding  curves  with  the  sum 
of  phase  angles  0^  and  0^  equal  to  180  degrees  are  the  points  of  the  locus 
required.  The  curve  drawn  through  these  points  is  the  required  locus.  It 
is  the  required  locus  because  each  point  on  it  satisfies  equation  4.16. 

The  determination  of  the  gain  can  be  carried  out  then  by  the  criterion  in 
the  root- locus  method.  The  magnitude  equation  is 


K  "  |CH<r  +  j of ) |  e"Ts  (4.17) 

which  will  give  the  gain  K  of  any  point  on  the  required  locus. 

The  required  roots  are  the  points  where  the  value  of  gain  K  is  equal 
to  one.  The  gain  limit,  which  is  the  maximum  gain  below  which  the  system 
can  have  damped  oscillation,  can  be  found  readily.  This  is  the  value  of 
gain  at  the  point  of  intersection  of  the  required  locus  and  the  imaginary 
axis.  This  value  is  found  to  be  1.21. 

This  given  function  transfer  has  an  infinite  number  of  roots,  which 
can  be  obtained  by  constructing  additional  branches  of  required  locus. 

In  general,  the  number  of  influential  roots  depends  on  the  relative 
magnitudes  of  the  dead  time  T  and  the  other  time  constant  in  G'(s). 

The  limit  cycle  occurs  at  CL/ -  1.67  rad/sec.  and  k  =  1.21  where 


k  =  kfk 


K1  = 


4Y 
n  X 


2_ 

JT  X 
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k  =  1.5 


K  =  1.21  = 
X  = 


n  X 
2  x  1.5 


«  x  1.21 


Kt  n  X 
=  .788 


x  1.5 


37 


5.  System  with  Dead  Time, 

5*1,  Transfer  function  associated  with  pure  dead  time 

If  the  input  to  component  with  a  pure  dead  time  characteristic  is  f(t) 
u(t),  then  the  output  of  the  component  is  f(t  -  T),  where  T 0  is  the 
dead  time/  The  transfer  function  associated  with  pure  dead  time  is,  there¬ 
fore, 


G(s)  - 


-TS  ,  N 
e  F(s) 

F(s) 

-TS 


For  sinusoidal  input  at  frequency  clt 
in  equation  5.2  becomes 

i 

r» 


(5.1) 

(5.2) 

the  transfer  function  given 


G(j  ur)  =  (5.3) 

The  gain  is  thus  seen  to  be  independent  of  the  frequency  and  equal  to  unity, 
and  the  phase  shift  is  seen  to  be  directly  proportional  to  frequency. 

5.2  Frequency  response  of  a  system  with  dead  time. 

It  is  of  interest  to  compare  the  frequency  characteristics  of  a  unit 
with  pure  dead  time  T  and  a  unit  with  a  first-order  lag  transfer  function 
with  time  constant  equal  to  T.  That  is, 

-TC 

Gx(s)  =  e  (5.4) 


and 

c2<s>  =  tsVt  <5'5) 

The  gain  associated  with  G^CjcoO  and  GjCj^cO  are 

|Gl(j^)|  =  1  (5.6) 

r  ■  5n"j 

j  G2(j  os  )|  =  V 1  +  cal  (5.7) 


and  their  respective  phase  angles  are 


0l  =  -a/  T 


(5.8) 
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Fig.  5-1  Bode  diagram  for  a  pare  dead-cine  unit  and  a  ?i. rec¬ 
order  lag  unit  with  time  cor.r.t,  i  T  :d:  Vtto 
equal  to  the  dead  time  of  tr 
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and 


0?  =  -  arc  tan OJT  (5.9) 

The  Bode  diagram  for  the  two  units  are  plotted  as  shown  in  Fig.  5.1,  where 
it  can  be  seen  that  for  a?T  1,  the  units  have  nearly  the  same  effect  on 
a  given  input  signal.  Generally,  however,  since  for  higher  frequencies  the 
dead-time  unit  exhibits  more  phase  lag  and  more  gain  than  the  first-order 
lag  unit,  the  introduction  of  the  dead-time  unit  in  cascade  within  a  given 
closed-loop  system  has  more  adverse  an  effect  on  the  stability  of  the  system 
than  does  the  introduction  of  the  first-order  unit. 

Dead  time  has  the  characteristic  of  a  pure  phase  shift  since  the  magni¬ 
tude  is  always  one.  Seen  easily  from  the  graphs  drawn  in  Fig.  5.2  that  no 
other  dynamic  element  builds  up  phase  lag  so  much  as  dead  time.  In  combina¬ 
tion  with  other  dynamic  elements  the  phase  lag  of  dead  time  is  simply  added 
to  obtain  the  total  phase  lag. 


Fig.  5-2.  Dead-time  element. 

A  useful  approximation  for  dead  time  T  is  the  following  function  from 
the  Pade  list, 


,fos  _  2  -  TS 
G(S)  -  2  +  TS 


(5.10) 
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The  magnitude  of  this  function  is  everywhere  1.0  and  the  phase  lag  cor¬ 
responds  very  closely  to  dead  time  up  to  frequencies  given  by  a/I  =  0.5. 

In  addition,  this  function  may  be  easier  to  incorporate  into  complex  system 
analysis . 

The  destabilizing  effect  of  the  transportation  lag,  as  dead  time  is 
also  known,  and  clearly  illustrated  in  a  Nyquist  diagram.  For  example,  if 
a  unit  with  a  pure  dead-time  characteristic  given  by 

Gj (S)  =  e'TS  (5.11) 

is  connected  in  cascade  with  a  component  with  transfer  function 

ys>  -  s(s2+n  (5-i2> 

as  in  Fig.  5.3,  the  Nyquist  diagram  for  the  system  for  various  values  of 
the  dead-time  appears  as  in  Fig.  5.4.  As  can  be  seen  in  this  graph,  for 
small  values  of  transportation  lag  the  introduction  of  the  dead-time  unit 
merely  results  in  a  decrease  in  relative  stability,  but  for  large  values  of 
transportation  lag  the  system  can  be  unstable. 

The  effect  upon  stability  of  the  introduction  of  a  transportation  lag 
in  cascade  within  a  closed-loop  system  of  the  type  illustrated  in  Fig.  5-5 
(a)  can  be  determined  from  the  Nyquist  curve  for  the  original  system  in  the 
following  manner.  The  characteristic  equation  of  the  original  system  is 

1  +  Ga(8)Gb(s)  =  °  (5.13) 

and  that  of  the  system  which  is  obtained  by  introducing  a  transportation 

lag  T  within  either  the  forward  path  or  the  feedback  path  is 

1  +  Ga(s)Gb(s)e"TS  =  0  (5.14) 

Equation  (5.14)  can  be  rewritten  as 

Ga(s)Gb(s)  =  -e+Ts  (5.15) 

or,  in  terms  of  frequency  response,  as 

Ga(j^)Gb(j^)  =*ej"/T  (5.16) 
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R  + 


®r-~ 


Fig*  5~3o  Block  diagram  of  a  simple  feedback 
system  with  transportation  lag. 


Fig,  5-4,  Nyquist  diagram  for  the  system 
illustrated  in  Fig*  5-3, 
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(a)  (b) 


Fig.  5-5  Block  diagram  of  a  single  loop  feedback  system 


Fig.  5-6  Stability  locus  in  the  G  G  -plane  for  a 

Q.  D 

system  of  the  type  shown  in  Fig.  5-5(b). 


Fig.  5-7  Determination  of  phase  margin  in  a 
system  with  transportation  lag. 
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The  locus  of  G  (jcc/')G,  ( j  ckT  )  is  the  Nyquist  curve  for  the  original  system; 

a  b 

equation  (5-16)  states  that  the  effect  of  introducing  transportation  lag  T 

is  to  shift  the  critical  stability  point  in  the  G  G,  -plane  from  -1  +  jO  to 

a  d 

-coso/T  -  jsimxrT.  That  is,  the  critical  point  becomes  a  critical  locus, 

in  the  form  of  the  unit  cycle  in  the  G^G^  -plane,  as  shown  in  Fig.  5-6. 

The  characteristic  equation  (5-16)  is  satisfied  only  if  the  locus  of 

G  ( jc</)G,  ( j  cu? )  intersects  the  stability  locus  in  such  a  manner  that  the 
a  b 

frequency  on  the  Nyquist  curve  and  the  frequency  on  the  stability  locus  are 
the  same  at  the  point  of  intersection.  If  the  two  frequencies  are  not  equal 
at  the  point  of  intersection,  a  measure  of  relative  stability  is  afforded  by 
the  stability  locus.  Thus,  if  the  Nyquist  diagram  for  a  given  system  has  the 
appearance  shown  in  Fig.  5-7,  the  phase  margin  can  be  determined  as  follows: 

1.  Determine  the  cross-over  frequency  ,  at  whichthe  Nyquist 

curve  crosses  the  unit  circle. 

2.  Locate  the  point  on  the  stability  locus  for  ccr  =  clS^  . 

3.  Draw  a  radial  line  segment  from  the  origin  through  each  of 
these  points. 

4.  Measure  the  angle  between  these  line  segments.  This  angle  is 
the  phase  margin,  with  the  same  positive  sense  as  in  the  case  of  a  system 
without  transportation  lag. 
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6.  Transportation  lag. 


6.1.  General 

The  phenomenon  of  dead  time,  sometimes  called  transportation  lag  or 
time  delay  or  ideal  time  lag,  can  best  be  categorized  as  a  lack  of  change 
in  output  for  a  specified  interval  of  time  following  a  change  in  input, 
this  is  illustrated  in  Fig.  6.1. 


Fig.  6.1 

Transportation  lag  may  be  encountered  in  process  involving  pressure, 
temperature,  and  also  may  be  found  in  process  involving  chemical  reactions 
when  a  finite  time  must  elapse  before  a  reaction  begins  to  occur. 

Occasionally  in  analyzing  engineering  systems,  it  becomes  necessary  to 
simulate  time  delay  (or  transportation  lag).  For  example  in  studying  trans¬ 
mission  processes  through  hydraulic  or  pneumatic  lines  or  in  studying  the 
motion  of  radio  active  particles  in  neuclear  reactors  it  is  necessary  to 
produce  an  output  voltage  which  is  identical  with  an  input  voltage  but  dis¬ 
placed  in  time  by  a  specified  amount. 

Transportation  lag  may  be  expressed  by  the  following  mathematical  re¬ 
lationship:  Given  f(t)  for  t  ^rO,  the  function  f(t  -  T)  differs  from  the 

function  f(t)  by  the  dead  time  T,  where  f(t  -  T)  is  defined  to  be  zero  for 
t  <  T. 
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A  direct  and  accurate  method  for  generating  such  time  delay  is  to 
employ  a  magnetic  tape  recorder  in  which  the  space  between  the  recording 
and  reading  heads  is  adjustable,  by  suitably  controlling  the  relative  loca¬ 
tion  of  the  two  heads  and  the  speed  of  the  tape  a  wide  range  of  time  delay 
can  be  simulated. 

Approximate  methods  which  do  not  require  special  auxiliary  equipment, 
once  an  analog  computer  is  available,  are  based  upon  the  Laplace-transforma- 
tion  relationship 

/[ / (t  -  Tj]  =  Ffs j e 5T=  U[ f«>] e'r  (6.1) 


-  sT 

since  e  cannot  be  expressed  as  a  rational  function,  the  transfer  func¬ 
tion  for  dead  time  can  be  simulated  only  approximately,  thus  the  remaining 

c 

methods  depend  on  an  approximate  series  or  technique  of  network  synthesis,"' 
The  Taylor  series  expansion  of  this  function  is 


/  -  TS  -t 


(TSX 

2  I 


(TS)d  ,  (Tsf 
3  1  '  4  / 


f 


(Tjf 

a  / 


(6.2) 


The  generation  of  the  terms  of  this  series  would  require  successive  differ¬ 
entiation  with  attendent  noise  problems.  The  rate  of  convergence  of  equa¬ 
tion  6,2  is  further  relatively  slow,  particularly  for  large  values  of  T, 
or  high  frequencies.  For  this  reason  the  Taylor  series  is  not  well  suited 
to  the  generation  of  transportation  lag  for  problems  involving  high  fre¬ 
quencies  or  long  time  constants. 

6.2  Pade  approximation 

A  more  practical  approach  to  the  generation  of  the  transportation  lag 

9 

function  is  to  employ  the  so  called  Pade1  approximation.  From  equation 
6.1  the  transfer  function  of  transportation  lag  is 

(C =fST _  'A  - ^2. -  (6.3) 
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substituting  s  =  jar  to  find  the  frequency  response  gives 


eJUJT=  Cos«jT  -  jS.nco'T 

E c 


thus 


ft 


=  Jc&S  2 as  T  -f  S>'n*co  / 


■=  fan  '(- 


5  is?  oo_  7~ 

Cos  as  T 


■) 


■  / 

—  co  ~ 


(6.4) 

(6.5) 


from  the  above  equations  the  frequency  transfer  function  of  the  trans¬ 
portation  lag  is  characterized  by  an  amplitude  ratio  which  is  unity  for  all 
frequencies  and  a  phase  angle  which  varies  linearly  with  frequency  over  the 
desired  range  of  interest.  From  equations  6.4  and  6.5  it  is  seen  that  the 
quotient  of  any  pair  of  complex-conjugate  polynominals  has  an  amplitude  of 
unity  and  also  possess  considerable  phase  shift. 

In  the  Pade'  approximation  the  exponential  is  approximated  as  the  ratio 
of  two  polynominals  in  s,  according  to 


 P„jS) 


(6.6) 


where  subscripts  m  and  n  refer  to  highest  power  of  s  in  the  numerator  and 
denominator,  respectively.  Expressed  in  its  basic  form,  the  Pade*  approxi¬ 
mation  is 


<S?X -  Li'm  (2Ll 

(o+v)  —  oo  GrutV  (x.) 


where 


l/X  I s(v-‘)X* 

(O-H/)1.  (.u  -tv)(u-n/->)2  / 


V(v-_ 0(w-2)  -  -  •  .  •  /  x  ^ 

(<-/+  I/)  (o+  is-/)  (CJ  +  I )  1/  I 


(6.7) 


(6.8) 
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-4-  *  « 


£•</,(/ (V —  / 


ux 

(Vtu)! 


a  (a- / )  x.1 


(v*  u)  (  (T-, <■  u-i )  2  I 


r  ,\U  iK2z  81  •  -  ■  £ '  /  *  * 

■+  (  ')  JpTuJcyTcT^r)  •  •  -  •  ■  (^-t- 1)  u  j 


(6.9) 


the  convergence  for  this  series  expansion  is  quite  rapid  and  often  values 
of  u  and  v  of  2  will  give  good  accuracy  for  short  time  lag  and  low  frequenc¬ 
ies. 


-  sT 

The  Pade'  approximation  of  e  for  u  =  v  =  2  and  if  x  =  -sT  is 


-Sr 

/  -  §-(sr)  T~r,  (s-rr 

(6.10) 

—  ' — v 

/  t  $  CST)  +  -3_  (ST)Z 

e  ■ ST  ^ 

S2Tz  -  657  ■+  )  2 

(6.11) 

S*T2  +  6  ST  -t  1 2 

the  higher  the  powers  of  the  numerator  and  denominator,  the  more  terms  of  the 

series  can  be  matched  and  the  smaller  the  error  in  the  time  delay  function. 

Equation  6,11  is  the  second  order  Pade1  approximation.  The  fourth  order  Pade* 

-  sT 

approximation  (u  =  v  =  4)  for  e  is 

-sr  (ST)*  ~20(sr)3  +  1 80(*7J  ~  840(t>l)  +  /C80 

&  ~  Z  - ; - - - 5 - (b .  12 ) 

{ST)  4  2o(sr)3  -r  l So  (.st)3-/-  S40(ST)  -t  1680 

Higher  order  pade'  approximations  provide  improved  accuracy  and  fre¬ 
quency  response.  It  is  noted  that  the  order  of  the  characteristic  equation 

-  sT 

is  raised  by  the  order  of  the  approximant  used  to  simulate  e  .  If  the 
second  order  approximant  is  used,  two  additional  roots  are  introduced  into 
the  characteristic  equation. 

6.3  Simulation  of  transportation  lag 
Equation  5.12  can  be  rewritten  as 


_  20_  /_80  840  ,  I08o 

J[-f(t-TJ]  =  7  '  sr  (srf  (srp  _  cjrr  (6vl3) 

7  ST  '  C*V*  ( 5T )3  r  (STj* 

since  the  differential  operator,  D,  and  Laplace  operator  are  interchangeable 


'‘8 


when  initial  conditions  are  zero,  equation  6.13  can  be  expressed  as 


/  (*  -  T) 
-f(t) 


2_Q  l_8o  _  &4Q  +  t±£Q 
'  D~f  '  (.o  ~)2  (o  r)*  (or)* 


(6.14) 


/  +  20  +  •8°  8±P  +  1680 

'  D  T  (or)2  (OTp  (Dr)*? 


rearrange  equation  6.14  and  we  get 


i«-T)=  -f(t)  -  +  A^f(t-r)}-^ 


-  f  (t  -  rjJ  -O. 


(6.15) 


according  to  equation  6.15,  the  analog  computer  circuit  can  be  simulated 
as  shown  in  Fig.  6.2, 

The  time  delay  T  in  second  is  determined  by  the  magnitude  of  the  input 
resistors  to  the  four  integrators.  Such  a  circuit  works  well  only  in  a 
limited  frequency  range  (direct  current  to  approximately  50  cps).  In¬ 
stability  results  if  the  variation  of  the  input  signal  f(t)  is  too  rapid. 

It  is  therefore  not  feasible  to  apply  a  step  function  to  this  circuit.  In 
case  f(t)  is  not  equal  to  zero  at  the  commencement  of  the  computer  run,  it 
is  necessary  to  apply  appropriate  initial  conditions  to  the  integrators. 

The  transportation  lag  devices  described  thus  far  are  limited  to  applica¬ 
tions  in  which  the  delay  time  T  remains  constant  throughout  the  computer 
run. 

The  method  of  simulating  transportation  lag  by  the  Padef  approxima¬ 
tion  yields  what  may  be  termed  optimally  flat  phase  characteristics  at  the 
origin  but  does  not  fit  the  phase  characteristics  as  well  over  an  extended 
frequency  range.  The  major  drawback  to  the  use  of  the  Pade1  approximation 
to  simulate  transportation  lag  is  their  response  to  a  step  input  or  other 
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Pade  approximac ion. 


input  having  high  frequency  components,  since  the  amplitude  ratio  is 
always  unity  but  the  phase  shift  is  approximated  only  for  low  C(f T  pro¬ 
ducts.  The  phase  shift  characteristics  of  certain  of  the  Pade'  approximation 
are  shown  in  Fig.  6.3.  As  a  result  of  the  inability  of  this  approximation  to 
linearly  phase  shift  high  frequencies,  the  output  of  an  approximant  rings 
badly  for  high-frequency  inputs.  This  behavior  is  illustrated  in  Fig.  6.4. 

The  actual  location  of  the  transportation  lag  transfer  function  in  the 
forward  path  is  immaterial  and  it  may  be  placed  at  the  end  of  the  path  in 
the  computer  where  the  high  frequency  components  have  been  attenuated  to  the 
greatest  possible  extent. 

In  the  great  majority  of  cases  the  fourth  order  Pade'  approximant  is 

more  than  adequate  for  the  simulation  of  transportation  lag.  In  cases 

where  linear  phase  shift  of  higher  ccrT  products  than  those  shown  in  Fig. 

6.5  is  desirable,  higher  order  Pade1  approximants  or  more  complex  curve 

-  sT 

fitting  techniques  may  be  used  to  simulate  e  .  The  phase  shift  character¬ 
istics  of  the  sixth  order  and  eighth  order  Pade1  approximants  are  shown  in 
Fig.  6.5.  On  the  basis  of  this  result  it  may  be  inferred  that  the  order 
of  the  Pade'  approximation  may  be  increased  to  increase  the  magnitude  of 
the  cw  T  product  that  is  linearly  phase  shifted  without  distorting  the 
phase  shift  of  low  t  products.  Because  the  Pade1  approximation  may  be 
expressed  in  a  general  form  from  which  any  order  approximant  may  be  develop¬ 
ed,  it  is  felt  that  the  use  of  curve  fitting  techniques  requires  a  needless 
expenditure  of  effort  without  any  tangible  benefit. 

6.4  Generation  of  time  delay  from  magnetic  tape  recorder. 

The  transportation  lag  generated  by  using  the  magnetic  tape  recorder 
is  very  accurate,  and  it  saves  a  lot  of  time  and  labor  for  simulation  by 
using  amplifiers. 
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WT  Product  (radians) 


Fig.  6-3  Phase  Shift  Characteristics  of  the  Fade 

Approximat ion 
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(a)  delay  time  ^  1»0  second 


(b)  delay  time  ^  10  seconds 

Fig.  6-4  Curves  obtained  from  Pade  approximation 

l 
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Phase  shift  (radian) 


Fig.  6-5  Phase  shif  Oharacterist ics  of  the 
Pade  Approximation 
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The  magnetic  tape  recorder  used  has  the  following  specifications: 
Ampex  Model  FR1107  magnetic  tape  recorder 
Tape  speed:  3-3/4,  7-1/2,  15,  30,  and  60  inches  per  second. 

Input  impedance:  100,000  ohms  unbalanced  to  ground . 

Input  level:  1  volt  rms  recommended  for  normal  recording  level. 

Output  impedance:  125  ohms  75  /yf,  unbalanced  to  ground. 

Output  level:  1  volt  rms  normal,  across  10,000  ohms  load  impedance  at 

normal  recording  level. 

The  response  curves  obtained  from  the  magnetic  tape  recorder  are  as 
shown  in  Fig.  6.6.  Compare  this  with  that  obtained  from  analog  computer 
as  shown  in  Fig.  6.4,  it  is  obvious  that  the  magnetic  tape  recorder  is 
most  desirable  for  simulating  the  time  delay,  since  it  gives  a  more  ac¬ 
curate  result. 
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(a)  delay  time  ^  .4  second 


(b)  delay  time  ^  1.7  seconds 
Fig. 6-6  Curves  obtained  from  magnetic  tape  reoorder 

I 
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7 .  System  simulation. 


7.1.  Analog  computer  simulation 

The  block  diagram  of  the  system  to  be  investigated  is  as  shown  in 

Fig.  7.1. 


Fig.  7.1.  The  block  diagram  of  the  relay  servo  system. 

The  simulation  of  the  transportation  lag  is  shown  in  Fig.  6.2.  The  simula¬ 
tion  of  the  ideal  relay  is  as  shown  in  Fig.  7.2. 


Fig.  7.2.  The  simulation  of  the  ideal  relay. 

From  the  block  diagram  of  the  system,  the  analog  computer  circuit 

diagram  can  be  simulated  as  shown  in  Fig.  7.3. 

Because  of  its  convenience  and  accuracy,  the  magnetic  tape  recorder  is 

_  s  t 

selected  for  the  generation  of  transportation  lag  (e  )  through  the  experi¬ 


ments. 
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Fig.  7-3  Operational  Diagram  of  Analog  Computer  setup  for  the 
block  diagram  of  the  relay  servo  system  of  Fig.  7-1. 


7.2.  Digital  computer  program 


The  block  diagram  of  the  system  is  shown  below. 


f(t)  +  _  +/-N  5(t) 


E  (IT) 


Delay 

Relay 

b(t)  ^ 

IDl 

EDELAY^ 

V 

B(IT) 

CTACH 


DA  DELAY 


ADELAY 


DELAY 

ID2 


Fig.  7.4.  The  block  diagram  of  relay  servo  system. 

The  variable  names  used  in  the  program,  E(IT),  EDELAY,  etc.,  are  shown 
below  their  mathematical  counterparts.  Note  that  the  fixed  point  IT  or  the 
floating  point  T  indicates  the  number  of  increments, starting  with  IT  =  1 : 


zero  time 

i 


r 


.time  =  DELAT  time  =  2  DELAT,  etc. 


IT  =  1 
T  =  1. 
E  (1 ) 


IT  =  2 
T  =  2. 
E  (2) 


etc.  etc.  etc. 

The  method  of  computation  is  simply  to  evaluate  the  analytic  expression  for 
the  response  at  each  increment  point.  The  response  is  always  supposed  to 
have  started  at  the  time  of  the  last  relay  operation  (TF  in  increment  measure) 
with  initial  condition  equal  to  the  final  values  of  A(IT)  =  AF,  and  DA(IT)/ 
CTACH  =  DAF  at  the  end  of  the  previous  response  segment: 
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new  response  computation  starts  here 
based  on  relay  operation  at  IT  =  TF 


J 


^if  computation  at  this  point  shows  change  of  sign  of  EDELAY  from  previous 
lue,  the  relay  operates  here,  TF,  AF,  and  DAF  are  recorded  as  the  current 


va 


value  of  IT,  A(IT)  and  DA (IT) . 
Analytic  response. 

The  plant  transfer  function  is 

A(s)  / 


B(V 


5(3  t  <*) 


e(s)—  S2/l(s)  -t  o(SA<s) 

and  the  differential  equation  describing  the  plant  is 

b(t)=  +  <* 

with  initial  conditions  a(0)  and  a(0)  (AF  and  DAF  in  program). 
The  input  is  a  step  V 

^  =  5a/l(s)  -  Sa(o)  -  d(o)  +  oCSACs)  -  <*cl(o) 

A  _  Vt  Sjq(Q)  ±JLL$LiE2 

'•  S2(S-h  cx) 

This  can  be  expressed  into  partial  fractions 

A(s)  =  JL.  +  .JL  _  X(s  ■<■*)  +  vs(s-f«j+2s2 

s 


5  S  +<* 


S2  (S  ■+  ot  ) 
2 


(7.1) 

(7.2) 

(7.3) 


(7.4) 


(7.5) 


equating  the 

coefficient  of  s 

5^  CS-+  * ) 

S2  : 

V  -+■  z  = 

a.  (o) 

(7.6) 

5'  : 

X  ■+  Yoi 

d{o)  -h  o(CUo) 

(7*7) 

S°  : 

x°<  — 

V 

(7.8) 
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X  =2.  ^  -  =  V/ALPHA  in  program 

^  I  (2(0)  - {-  o(.<Z(o)  -  )()  =  (DAF  +  ALPHA-X) /ALPHA 

^  =  aco)  "  Y  =  AF  -  Y 

from  equation  7.5,  the  time  solution  is 

-at 

act)- xt  +  y  -tzG 

=  X*TAU  4  Y  4  Z*EXPF ( -ALPHA *TAU)  in  program 

a(±)-  x 

=  X  -  Z*ALPHA*EXPF(~ALPAH*TAU)  in  program 
t-^™^TAU  is  time  from  last  relay  operation 
Program  accuracy. 

1) .  The  time  delays  can  be  introduced  only  as  an  integral  number 
of  increments  (each  of  magnitude  DELTAT) . 

2) .  The  only  actual  inaccuracy  is  the  fact  that  the  relay  can  only 
operate  at  integral  values  of  IT*.  The  accuracy,  for  a  given  increment  size, 
DELTAT  could  be  fairly  easily  improved  by  introducing  an  interpolation  pro¬ 
cedure  between  computation  points.  The  same  technique  could  be  used  to 
introduce  non-integral  number  of  the  delays. 

The  plant  dynamics  can  be  easily  changed  by  re-working  the  partial 

fraction  expression  and  the  expression  for  a(t)  and  a(t). 

The  program  for  the  relay  servo  system  is  shown  below.  And  one  case 

of  its  output  (t,  =  .8  second,  k  =  2,0)  is  attached  to  this  program, 
d  t 


*This  error  may  be  regarded  as  a  variability  in  the  forward  path  time  delay, 
so  that  the  actual  delay  may  lie  between  the  desired  ID1  and  ID1+1. 
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.JOB  CHI ANG  ,  MAX.  TIME  3  MIN. 

PROGRAM  TLRSOJ 

RESPONSE  CALCULATIONS  FOR  RELAY  SERVO  SYSTEM  WITH  DELAYS  ID1.ID2  AND 
ID3  IN  THE  FORWARD, FEEDBACK  AND  TACHOMETER  PATHS.  RESPECTIVELY. 
SOLUTION  IS  OBTAINED  FOR  A  TOTAL  TIME  OF  TMAX  SECONDS  IN  INCREMENTS 
OF  DELTAT  SECONDS.  IN  THE  PROGRAM  BOTH  T  AND  IT  ARE  THE  NUMBER  OF 
TIME  INCREMENTS  PLUS  ONE.  THAT  IS,  THE  IN  I TI  ALr’  T IME  IS  T=1..  THE  THREE 
DELAYS  ARE  MEASURED  BY  THE  NUMBER  OF  TIME  INCREMENTS.  THE  SYSTEM 
FORCING  FUNCTION  IS  SPECIFIED  BY  PROGRAM  STATEMENT  NUMBER  1,  AND  TWO 
FOLLOWING  CAROS.  v 
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INPUT  DATA. 

IN 

101,102, 

V 

ALPHA 

CTACH 

TMAX 

OELTAT 

INCPRNT 

INCPRNT  IS  THE 
V  IS  THE  RELAY 
ALPHA  IS  THE 
CTACH  IS  THE 


ON  CARD 
ON  CARDS 
ON  CARD 
ON  CARO 
ON  CARD 
ON  CARD 
ON  CARO 
ON  CARO 
NUMBER  OF 
VOLTAGE 
PLANT  TIME 
TACHOMETER 


1 

2 


RIGHT  JUSTIFIED  iM  COLUMNS  1  AND  2 


THRU  4.  RIGHT 


COLUMNS 

5.  IN  COLUMNS  1  THRU  10  WITH  DECIMAL 

6.  IN  COLUMNS  1  THRU  10  WITH  DECIMAL 

7.  IN  COLUMNS  1  THRU  10  WITH  DECIMAL 

8.  IN  COLUMNS  1  THRU  10  WITH  DECIMAL 

9.  IN  COLUMNS  1  THRU  10  WITH  DECIMAL 

10.  RIGHT  JUSTIFIED  IN  COLUMNS  1  AND  : 

INCREMENTS  BETWEEN  PRINT  OUTS 


1  THRU  3 

POINT 

POINT 

POINT 

POINT 

POINT 


CONSTANT 

CONSTANT 


1 02, 1 03, V, ALPHA, CTACH, TMAX, CELT AT, INCPRNT 


THE  MAXIMUM  NUMBER  OF  INCREMENTS  IS  1000.  THE  NUMBER  OF  INCREMENTS 
MAY  NOT  EXCEEO  900  IF  A  GRAPH  IS  REQUIRED. 

IN  IS  THE  NUMBER  OF  CASES  TO  BE  SOLVED.  NOTE  THAT  THERE  MUST  BE  ONE 
DATA  OECK  OF  NINE  CAROS  FOR  EACH  CASE.  THE  FIRST  CARD  CONTAINING  IN 
IS  NOT  REPEATED. PUT  IN=1  IF  ONLY  ONE  CASE  IS  TO  BE  SOLVED. 

EACH  SET  OF  DATA  CAROS  MUST  CONTAIN  THE  REQUIRED  GRAPH  FORMAT  CAROS. 
REMOVE  CALL  GRAPH  STATEMENT  AND  GRAPH  FORMAT  CARDS  IF  NO  GRAPH  IS 
REQUIRED. 

DIMENSION  A{ 1001 > ,B( 1001 ) ,E( 1001 ) ,DA( 1001 ) ,DE( 1001 } ,ET ( 1001 ) , 

1  TIMES! 1001 ) ,FIN(1001 ) » 8S  I 1001) 

READ  150,  IN 
FORMAT  (12) 

00  34  K  =  1 , IN 
READ  100,101, 

PRINT  200 
PRINT  201 

PRINT  202, 101 ,102, 103, V, ALPHA, CTACH, TMAX, DELTAT 
FORMAT  ( 1 3/13/ I3/F 10.1 /FI 0.1 /FI  0.1 /F 10.1 /F10. 1/12) 

FORMAT  {//) 

FORMAT ( 3 1 H  RELAY  SERVO  SYSTEM  INPUT  OATA  //) 

FORMAT  t  8H  ID1  =13/  8H  102  =13/  8H  103  =I3/8H  V  =F7.3/ 

1  8H  ALPHA  =F7 .3/  8H  CTACH  =F7.3/  8H  TMAX  =F7.3,5H  SEC./ 

2  ’  8H  DELTAT*F7. 3, 5H  SEC.) 

PRINT  203 
PRINT  204 

FORMAT  { 6X42HCALCULATED  RESPONSE  OF  RELAY  SERVO  SYSTEM  ) 

FORMAT  (1H1) 

T  =  1  . 

TF=  1  . 

A  (  1 ) =0. 

DAI  1 ) =0 . 

AF=0  • 

D AF=0 • 

EDELAY=0. 

INC=TMAX/DELTAT 
NUM= INC-2 
F=  1  . 

DF=0. 

FORMAT  I6X30HF0RCING  FUNCTION  *  A  UNIT  STEP) 

FINI IT)=F 
IFIT-1.)  7,2,7 
PRINT  205 
PRINT  200 
PRINT  206 
PRINT  207 

FORMAT  ( 7X4HTI ME5X7HF0RC ING5X6HSYSTEM5X5HPLANT4X5HPLANT7X4HRATE 
1  7X4HRATE) 


150 


100 

200 

201 

202 


204 

203 


1 

205 


206 
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3 

4 

5 

6 


208 


207  FORMAT  ( 1 6X8HFUNCT I0N4XSHERR0R6X5H I NPUT4X6H0UtPUT5X6H0UTPUT 
1  6X5HERROR/ ) 

E ( 1 ) =F 
ET ( 1 ) =F 
I F ( E ( 1 ) ) 3. 5. 3 
I F ( 101)5,4.5 
B  ( 1  )  =SI  GNF (  V, E( 1  )  ) 

GO  TO  6 
B ( 1 ) *0. 

TIME*0. 

0E( 1 ) =DF 

PRINT  208,  TIME,F,ET( 1 ), B( 1 ) , A ( 1 ),DA( 1 ),DE( 1 ) 

FORMAT  (6XF6.3,3XF8.4,4XF8.4,3XF6.2,3XF8.4,3XF8.4,3XF8.4) 

DO  33  I T=2 » I NC 
T  =  I  T 

TIME=(T-1 . )*OELTAT 
TIMES! I T ) =TI ME 
EDELAY1=EDELAY 
J  =  1 

GO  TO  1 

7  TAU= ( T-TF ) *DELTAT 
TAU  IS  TIME  FROM  LAST 
TF  IS  THE  TIME  OF  THE 
X=B( IT-1 ) /ALPHA 
50  Y*(DAF+ALPHA*AF-X) /ALPHA 
Z=AF-Y 

AF  AND  DAF  ARE  THE  PLANT  OUTPUT  AND  RATE  OUTPUT  AT  TIME  OF  LAST  RELAY 
OPERATION.  X, Y  AND  Z  ARE  THE  COEFFICIENTS  IN  THE  PARTIAL  FRACTION 
EXPANSION. 

A( IT)=X*TAU+Y+Z*EXPF( —AL  PH A* TAU ) 

DA (IT)*  X-Z* ALPHA *EXPF (- ALPHA*TAU) 

I F ( 102)9.8,9 
ADELAY=A( I T ) 

GO  TO  12 

IFUD2-IT+1  )10. 10, 11 
ADELAY=A( IT-ID2) 

GO  TO  12 
ADEL AY=0. 

I F ( I D3 ) 14,13,14 
DADELAY=DA(IT) 

GO  TO  17 

14  IF(ID3-IT+1) 15,15,16 

15  DADELAY=DA( IT-ID3J 
GO  TO  17 
D ADELAY*0. 

E( IT )=F-ADELAY-DADELAY*CTACH 


RELAY  OPERATION  (OR  FROM  TIME  ZERO) 

LAST  RELAY  OPERATION  (INCREMENT  MEASUREMENT) 


8 

9 

10 

1 1 
12 
13 


16 

17 


51 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 
210 

29 

30 

331 

332 

333 


GO  TO  (51,31 ),J 
I F ( ID1 ) 19, 18,19 


EDELAY=E( IT) 

GO  TO  22 

I F ( ID1-IT+1 ) 20.20, 21 
EDELAY=E( IT-ID I ) 

GO  TO  22 
EDEL AY=0. 

IF ( E DEL AY  124,23,24 
B ( I T ) =0 . 

GO  TO  25 

B(IT)=SIGNF(V,EDELAY) 

IF ( XMODF ( ( IT-1 ),40*INCPRNT> ) 27,  26,27 
PRINT  203 

IF (XMODF ( ( IT-1 ) , 10*INCPRNT) ) 29,28,2 9 
PRINT  210 
FORMAT  (/) 

I F(8( I T ) —  B ( IT-1 ) ) 30, 31 *3  0 
IFlIDl-IT+l) 332,331 , 332 
S=0. 

GO  TO  333 

S=ABSF(EDELAY)/(ABSF(EDELAY)+ABSF(EDELAY1) ) 
TF=T-S 

AF*A(IT)-S*(A( IT )— A( I T— 1 ) ) 


I 
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DAF  =  DA(  IT)-S*(DAUT)-DA{  IT-1)) 

C  THE  LAST  SEVEN  STATEMENTS  SET  THE  INITIAL  VALUES  FOR  THE  NEXT  SEGMENT 
C  OF  THE  CALCULATION  IF  RELAY  OPERATION  HAS  BEEN  DETECTED  BY  STATEMENT 
C  NUMBER  29.  THE  TIME  OF  SWITCHING  HAS  BEEN  FOUND ‘BY  LINEAR 
C  INTERPOLATION  AS  HAVE  AF  AND  DAF.  F  SHOULD  BE  CONTINUOUS  EXCEPT  AT 
C  TIME  ZERO.  '* 

T IME=( T-S-l • ) *DELTAT 
PRINT  209. B( IT) .TIME 

209  FORMAT  (22H  THE  RELAY  SWITCHED  TO  F6.1.10H  AT  TIME  =  F7.3.5H  SEC.) 
C  WE  HAVE  NOW  TO  CORRECT  A ( IT) . 8( IT )  AND  EC  IT)  FOR  SWITCHING  AT  TF. 

T IME= (T- 1 • )*OELTAT 
J=2 

TAU= (T-TF)*DELTAT 
X  =  B ( IT ) /ALPHA 
GO  TO  50 

BSC IT)=B( IT)«0. 125/ABSFC  V)+0. 125 
ETC IT)=F-A( IT) 

DEC IT)=DF-DA(IT) 

IFCXMODFC ( IT-1 ) , INCPRNT) ) 33.320. 33 
PRINT  208, T I  ME, F, ETC  IT), BC IT ) , AC  I T ) ,0A( I T ) , DEC  I T ) 

CONTINUE 

( NUM.ET ,DE,8) 

( NUM, A  , DA , 8 ) 

(NUM.TIMES.A, 8) 

(NUM, TIMES, DA, 8) 

(NUM, TIMES, BS, 8) 

(NUM,  TIMES, ET, 8) 


31 


320 

33 


34 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

PRINT 

END 

END 


GRAPH 

GRAPH 

GRAPH 

GRAPH 

GRAPH 

GRAPH 

203 


4 

0 

0 

0 

.5 

1. 

.02 

30. 

0.05 

10 

0  +O.OOE+00  ♦O.OOE+OO 

02  EDOT  VS.  E 
02  CHI ANG 
0  ^  0 

0,  +0.00E+00  +0.00E+00 

02  COOT  VS.C 
02  CHI ANG 
0  0 

0  +4.00E+00  +5.00E-01 

02  C  AND  CDOT  VS. TIME 
02  CHI ANG 

1  0  C 

2  0  CDOT 


0 

0 


.5 

1. 

.02 

30. 

0.05 

10 

0  ♦0. OOE+OO  +0.00E+00 

02  EDOT  VS.  E 
02  CHI ANG 
0  0 

0  +0.00E+00  +0.00E+00 


00 

00 

08 

08 

1 

1 

02 

00 

00 

08 

08 

1 

1 

02 

04 

O 

o 

08 

08 

1 

1 

02 

00 

00 

08 

08 

1 

1 

02 

00 

00 

08 

08 

1 

1 

02 
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02  COOT  VS.  C 
02  CHI ANG 
0  0 

0  +4.00E+00 

02  C  AND  COOT  VS. 
02  CHI ANG 


+5.00E-01 

TIME 


04  00  08  08  1  1  02 


1 
2 
2 
3 

16 

0 

0 

.5 

1. 

2.0 
30. 
0.05 
10 
0 


c 

COOT 


+0.00E+00 
E 


02  EDOT  VS. 

02  CHI ANG 
0  0 

0  +0.00E+00 

02  CDOT  VS.  C 
02  CHI ANG 
0 

+4.00E+00 
C  AND  CDOT  VS. 


0 

0 

02 


♦0. 00E+00 

00 

00 

08 

08 

1 

1 

02 

+0.00E+00 

00 

00 

08 

08 

1 

1 

02 

♦5. 00E-01 

04 

00 

08 

08 

1 

1 

02 

02  CHI ANG 


TIME 


1  0  C 

2  0  CDOT 

2  0 

3  0  E 

20 

0 
0 
.5 
1. 

2.4 
30. 

0.05 
10 

0  +0.00E+00 

02  EOOT  VS.  E' 

02  CHI ANG 
0  0 

0  +0.00E+00 

02  CDOT  VS.C 
02  CHI ANG 
0  0 

0  +4.00E+00 

02  C  ANO  CDOT  VS. TIME 
02  CHI ANG 
IOC 

2  0  CDOT 

2  0 

3  0  E 


♦O.OOE+OO 

00 

00 

08 

08 

1 

1 

02 

♦O.OOE+OO 

00 

00 

08 

08 

1 

1 

02 

+5.00E-01 

04 

00 

08 

08 

1 

1 

02 

..END 
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RELAY  SERVO  SYSTEM  INPUT  DATA 


IDl  *  16 

ID2  *  0 

ID3  «  0 

V  ■ 

ALPHA  «  1. 

CTACH  »  2. 

TMAX  «  30. 
DELTAT* 


SEC. 

SEC. 


:f ; \ 

f  66  • 


500 

000 

000 

000 

050 


oo 

ON-JN 

cgmo 

o 

oo 

COJ- 

JtCOf- 

ro^rro 

coom 

cr'r- 

r-O'm 

r-co 

a: 

oo 

o— crco 

oro^ 

n- 

Lf>m 

-d*— • 

CMDCsJ 

rom-— 

mmm 

CNO 

*OCO 

csjm 

ujo 

oo 

Q'm-^o 

fOOO> 

csjro 

N-m 

ON-.^ 

ON-J* 

co^co 

r^-CNJr— 

ON 

♦-a: 

oo 

OCMfO^J* 

mo— 

o 

CNJ^O 

o— 

—  OCM 

cmocm 

O  —  CM 

o— 

CMOCM 

(MO 

<<* 

•  • 

•  tii 

#  •  • 

• 

•  ♦ 

•  • 

•  •  • 

♦  ♦  • 

•  ♦  • 

•  • 

•  •  ♦ 

t  • 

0CUJ 

1  1  II 

1  1 

1 

1  1 

i 

1  1 

1 

1  1 

i 

1  1 

1 

OO 

'Of'-j-r- 

eg  mo 

o 

OO 

oo«=r 

-3“C0N- 

m-=rm 

cocMm 

O'N- 

N-0'm 

f-CO 

K 

OO 

O— O'® 

om^f 

h- 

mm 

.=r— 

omcM 

mm— 

mmm 

ego 

OCO^J* 

cm  m 

UJO 

oo 

0'»r>^o 

moo' 

J- 

eg  m 

N-m 

ON^ 

ONJ1 

co^tco 

N-CM  — 

ON 

N-Ql 

oo 

otvjfojf 

mo— 

O 

com 

Or- 

f-OCsJ 

CMOCM 

O— CM 

o— 

CMOCM 

eoO 

<K 

a:3 

o 

I 

1 

1 

1  1 

l 

1 

1 

1 

UJ 

>»- 

<00 


oo 

.^mom 

mo— 

O' 

O'O' 

COO' 

m— cm 

—ooo 

mom 

O'LT) 

mNm 

mo 

►— 

oo 

O' COO— 

comm 

O 

CMJ* 

— 

O'J'N- 

m— O' 

eMesjm 

0'm 

mr^  cm 

—co 

i-o 

oo 

oo'mzr 

JtCNN- 

m 

om 

CMOCO 

cm  mo 

OCOO' 

com 

rnN^t 

N-O' 

za 

oo 

OOCMJt 

or-o 

o 

N-CO 

O' O' 

co  co  co 

OOO' 

COCO  O' 

oo 

as  co  o' 

OO 

<k 

-13 


ao 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

z 

00 

<✓> 

m 

m 

00 

oo 

00 

oo 

00 

oo 

UJ 

H 

o 

O' 

O 

CO 

N- 

CO 

o 

-3- 

o 

CO 

<✓> 

hh 

oooooooh-ooomo 

ooooooooo^oooo 

OOOOOONOOOOOO 

> 

ZD 

oocommmmoommmom 

mmommo'mmmmmmmN- 

mmm— mmja-mmmcomm 

m 

<a 

•  •  • 

♦  •  •  • 

•  •  ♦ 

•  t  «  • 

-IZ 

CM  1 

l  I  * 

m 

1  1  o 

CO 

1  1  1  O' 

c— 

1  1  CM 

m  i  i 

o 

Qlh 

— 

— 

*- 

cC 

HI 

u 

II 

ii 

ii 

11 

II 

II 

II 

II 

.  ii 

VOQ. 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Z 

UJC£ 


OO— O— 0*00  Hr- iD'O  — O' 
OOH-O'O^inKUOh-CslHO 


—  —  z  cm  —  z  n-  cr*  co  z  —  co  x  z 
h-  m  «  oo  cm  »-•  o  m  c\i  h  m  r-  o  . 
O'  m  h-m  o  h  n-  0s  r-  i—egm  ro  h- 


n-  n-  z  cmo  5:  n-  ro  n-  z  m  o 

N-  fs-  o  — <  CMV")  <-•  O  CM  N-  —  —  00 
m  —  O »-  C0m  h-OCMm  l-N-  O' 


-1 

HO 

oo  asa'N-m 

i  megm 

i  m 

eg— 

oo  — — 

p- 

ooo 

— — o  oo  o— 

o  oo 

UJh- 

ooa: 

•  •  • 

•h-  •  • 

OCH 

z 

>cc 

00  UJ 

— < 

< 

< 

< 

< 

< 

1  1  < 

<  1 

i  < 

<  1  1 

U-3 

m 

m 

m 

m 

m 

m 

m 

m 

m 

o 

< 

UJ 

00  II 

z 

‘  • 

• 

l 

• 

• 

l 

• 

• 

l 

• 

• 

i 

• 

• 

1 

Z  oo  oo  oooo  ooo  O 

oz  z*-«  ooooooooooooo 

0.0  *-»h-  OOhOOOOhOOOhO 

m—  oo  OO  OOOO  OOO  o 


oo  oo  ooo  ooo 

ooooooooooooo a 

OOHOOhOOOhOOOh 

OO  oo  ooo  ooo 


ooo  oo  ooo  oo 
ooooooooooooo 

OOOhOOhOOOhOO 

ooo  oo  ooo  oo 


UJh- 

acz 

•  *o 

•  •  •  «o  • 

•  «o  • 

•  «o  • 

•a  • 

•  *o 

•  ♦  «o 

•  •  •  o  • 

•  O  • 

•  •  Q  •  • 

a:o 

03 

—  UJ« 

—  —  UJ— < 

—  UJ— 

m*  —  UJ  —  < 

r-UJ^-i 

f— r—  UJ* 

— —  —  UJ 

—  —  —  LUr— 

»— LU^- 

r—  r—  UJr—  r— 

z 

U.U. 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

03 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

UJUL 

♦- 

h- 

►- 

h- 

H 

H- 

1— 

h- 

1- 

H 

H 

H 

>— « 

— « 

*—» 

M 

►— • 

<o 

3C 

X 

3: 

3: 

3 

s 

-IZ 

OOoOOOOOOOOOOoOO 

oomooooooomooom 

ooomoomooomoo 

3  — 

UJ 

OO  OOOO  ooo  o 

oo  oo  OOO  i 

ooo 

ooo  oo  ooo  oo 

OO 

Z 

om>omom>-omo>m 

om>om>omo>-mom>- 

omo>mo>mom>-om 

-la: 

H 

•  •< 
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7.3  Comparison  of  the  results  obtained  from  the  digital  computer  and  from 
the  analog  computer. 

The  result  obtained  from  the  digital  computer  for  the  system  with  time 

delay  .4  second  has  been  checked  with  analog  computer.  They  are  quite  in 

agreement  with  each  other  as  shown  in  Fig.  7.5  and  7.6. 

comparison  of  the  results 

by  digital  by  analog 

computer  computer 


overshoot  and  undershoot  3.5  2.6 

rise  time  5.9  5.7 

Because  the  digital  computer  can  give  an  accurate  and  exact  value  for 
any  desired  instant,  therefore,  the  digital  computer  is  selected  for  con¬ 
ducting  most  of  the  experiment  work,  and  the  analog  computer  is  used  to 
check  the  results  obtained  by  digital  computer.  They  are  quite  in  agreement. 
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Ob)  results  obtained  from  analog  computer 
Fig. 7-5  Output, velocity  of  the  system,  t^  -  . 


111 


4  second 


*7 1 


8.  Effect  of  transportation  lag  on  relay  servo  system. 


8.1.  General 

From  the  proceeding  sections*  we  know  that  the  time  delay  has  an  adverse 
effect  on  control  systems*  and  which  is  unavoidable  in  nearly  all  the  control 
systems.  Because  the  time  delay  may  basically  be  represented  by  an  infinite 
series  of  terms*  this  creates  an  infinite  number  of  roots. 

Later  on  in  this  section,  it  will  be  noted  that  the  time  lag  results  in 
successive  overshooting  and  undershooting  of  the  correct  final  value*  re~ 
presenting  serious  oscillations  about  the  final  position.  In  addition  to 
this  it  will  also  change  the  root  locus  and  phase  trajectories. 

8.2  The  frequency  response. 


From  Fig.  4.6(c)  we  have  the  transfer  function  with  removed 


G(s ) 


s  (s-i-  i ) 


-ST 


(8.1) 


where  K_  and  T  are  variables.  In  case  of  time  delay  equal  to  .5  second*  it 
is  found  that  the  tachometer  gain,  =  1.5,  will  compensate  the  system  to 
an  acceptable  optimum  operation,  this  will  be  seen’ in  section  9,  now  equation 
8.1  can  be  written  as 


G(s)= 

5(S-+‘) 


(8.2) 


The  Bode  diagram  for  the  above  transfer  function  is  shown  in  Fig.  8.1.  So 
far  as  the  linear  system  is  concerned*  the  gain  margin  for  .5  second  time 
delay  system  is  approximately  1.6  db,  but  for  the  undelayed  system*  it  is 
stable  for  all  gains  and  never  unstable. 

The  time  delay  has  introduced  additional  branches  into  the  diagram. 
These  branches  appear  due  to  the  periodic  nature  of  the  time  delay*  and  the 
increasing  phase  shift  introduced  at  higher  frequencies.  It  is  evident  that 
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Fig*  8-1  Bode  diagram  of  Eqc 


the  gain  margin  for  each  of  the  branches  increases*  This  is  analogous  to 
the  statement  that  the  roots  of  the  primary  branch  would  go  unstable  first 
if  the  gain  were  raised* 

The  gain  vs*  logos  curve  of  the  delayed  system  is  identical  to  that 
of  the  undelayed  system,  since  the  time  delay  introduced  no  additional  gain 
to  the  system*  The  phase  angle  vs*  log^u*  curve  of  the  delayed  system  is 
drastically  changed  since  the  equation  for  the  phase  angle  is 

-(90°  -h  OJ  -f  cot)  ~t~  tan' (<■  5co)  (8.3) 


The  CO T  term  in  the  above  equation  has  a  predominate  effect  to  the  phase 
angle,  it  introduces  additional  phase  shift  and  extra  branches* 

The  limit  cycle  occurs  at  OJ  =  1.7,  k  =  1*20* 

8,3  The  root  locus. 

The  root  locus  of  any  system  will  be  greatly  changed  when  a  transporta¬ 
tion  lag  was  introduced.  As  a  matter  of  fact,  for  an  undelayed  system  the 
transfer  function  in  Equation(8. 2)  becomes 

 A  5  (  s  '  &  y) 


C?  (  $) 


s  CS-hi) 


(8*4) 


no  limit  cycle  exists,  because  the  root  locus  never  intersects  the  imaginary 
axis  as  shown  in  Fig*  8.2,  which  shows  a  stable  system  for  all  gains* 


But  when  time  delay  is  introduced  this  situation  will  be  changed*  The  root 
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locus  for  Eq,  8,2  (,5  second  delayed  system)  is  shown  in  Fig.  8,3,  The 
construction  details  for  this  root  locus  were  described  in  section  4,4,2. 

Fig,  8,3  illustrates  the  effect  of  time  delay  on  the  root  locus.  It 
is  apparent  that  this  system  will  go  unstable  with  only  a  small  increase 
in  gain.  The  primary  branch  stability  limit  before  the  root  locus  corsses 
the  imaginary  axis  are  k  =  1,2, ,  and  OJ  -  1,76  rad/sec, ,  where  the  limit 
cycle  occurs,  this  almost  checks  with  Bode  diagram  where  k  =  1,20,  <xr  =  1,7 
rad/sec. 

Because  time  delay  introduces  periodically  phase  shift  of  360  degrees, 
so  there  presents  additional  (upper)  branches  when  at  high  frequencies.  It 
is  noted  that  for  the  same  gain  the  additional  roots  on  the  upper  branches 
are  attenuated  more  than  the  roots  on  the  primary  branch.  This  substantiates 
the  Bode  diagram  findings,  which  shows  that  the  roots  on  the  primary  branch 
phase  shift  would  cause  instability, 

8.4  Transient  response. 

Because  of  time  delay  introduced  in  the  system,  there  is  no  graphical 

method  that  can  be  easily  used  to  produce  the  transient  response  of  any 

system.  The  transient  response  of  the  delayed  relay  servo  system  (second 

order)  for  various  delay  time  and  various  k  (tachometei  gain)  are  shown 

t 

in  Fig.  8.4.  In  Fig.  8.4(a)  there  are  six  curves  under  the  same  taebo- 
meter  gain  (kfc  -  .02);  curve  (1)  is  an  undelayed  system,  the  rest  of  the 
five  curves  are  with  ,1,  .25,  ,5,  ,8  and  1,0  second  delay  time  respectively. 
For  1  second  time  delay  system  it  presents  an  overshoot  and  undershoot  of 
approximately  60  percent. 

Table  8,1  is  the  data  obtained  for  M  ,  t  ,  t  under  different  values 

Pt  P  r 

of  and  t^,  Plotting  these  data  into  curves  it  is  found  that  they  appear 

as  a  straight  line.  In  Fig.  8.5  are  shown  M  and  t  versus  delay  time,  in 

Pt  P 
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Fig.  8-3  Root  locus  of  relay  servo  system  with  time 

delay  of  .5  second. 
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-4  Cl,  Response  of  relay  servo  system  t,zO,  •!»  .25,  .5,  .8,  1.0  second  respectively. 
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Fig.  8 -4b  ,  Response  of  relay  servo  system  with  td=  0,  .1,  .25,  .5,  .8,  1.0  second 

respectively, k. -.8 


Table  8,1.  Data  for  first  overshoot  t  and  M  under  various  time  delays 

p  pt 

for  a  step  input. 


Kt 

Td 

M 

Pt 

t 

P 

t 

r 

amplitude  of 
limit  cycle 

.02 

.0 

1.1285 

3.50 

2.828 

no  limit  cycle 

If 

.1 

1.1870 

3.70 

2.95 

.08 

II 

.25 

1.2510 

4.10 

3.095 

.1625 

II 

.5 

1.3750 

4.60 

3.350 

.2575 

II 

.8 

1.5300 

5.10 

3.650 

.4910 

II 

1.0 

1.6500 

5.50 

3.850 

.5900 

.2 

.0 

1.0370 

3.35 

2.900 

no  limit  cycle 

II 

.1 

1.0870 

3.51 

2.950 

negligible 

If 

.25 

1.1630 

3.81 

3.080 

.057 

It 

.5 

1.2880 

4.30 

3.350 

.213 

II 

.8 

1.4370 

4.95 

3.650 

.388 

II 

1.0 

1.5310 

5.40 

3.850 

.49 

.4 

.0 

/ 

/ 

/ 

no  limit  cycle 

II 

.1 

/ 

/ 

/ 

II 

II 

.25 

1.0600 

3.70 

3.200 

.024 

II 

.5 

1.1650 

4.20 

3.420 

.125 

II 

.8 

1.3280 

4.80 

3.700 

.302 

II 

1.0 

1.4300 

5.20 

3.900 

.390 

.8 

.0 

/ 

/ 

/ 

no  limit  cycle 

II 

.1 

/ 

/ 

/ 

II 

II 

.25 

/ 

/ 

/ 

.015 

II 

.5 

1.0000 

3.70 

3.60 

.061 

If 

.8 

1. 1500 

4.39 

3.65 

.153 

If 

1.0 

1.2495 

4.71 

3.85 

.2495 

1.0 

.0 

/ 

/ 

/ 

no  limit  cycle 

If 

.1 

/ 

/ 

/ 

If 

II 

.25 

/ 

/ 

/ 

.014 

It 

.5 

/ 

/ 

/ 

.055 

II 

.8 

1.0550 

4.10 

3.70 

.137 

It 

1.0 

1.1550 

4.60 

3.85 

.205 
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t  (time  for  first  overshoot) 
P 


-  -  -  -  M  (magnitude  for  first 


overshoo  t.) 
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Fig,  8,6  is  shown  the  rise  time  versus  delay  time,  there  it  shows  that  the 
effect  of  time  delay  on  the  transient  response  is  to  increase  the  over¬ 
shoot  and  undershoot  and  cause  continuously  oscillations  and  make  the  rise 
time  longer,  it  turns  out  that  the  response  of  the  system  can  never  be 
settled. 

It  is  noted  that  M  ,  t  ,  t  increase  as  the  delay  time  increases, 

Pt  P  r 

But  when  we  increase  the  value  of  k^_  they  decrease.  And  to  a  certain  value 
of  k  some  of  them  become  acceptable  optimum  operation  this  will  be  seen  in 
section  9.2. 

In  Fig,  8.7  are  shown  the  plant  output,  rate  output  and  system  error. 
Fig.  8.7(a)  is  the  response  of  the  undelayed  system,  the  plant  output  set¬ 
tled  at  7  seconds;  Fig.  8.7(b)  shows  the  plant  output,  rate  output  and 
system  error  for  a  system  with  delay  time  of  ,5  second  and  k^  =  ,02,  the 
output  continuously  oscillates  with  an  overshoot  and  undershoot  of  approxi¬ 
mately  30  percent;  Fig,  8.7(c)  shows  the  plant  output  rate  output  and  system 
error  for  a  system  with  time  delay  of  one  second  and  k^  -  .02,  The  plant 
output  continuously  oscillates  with  an  overshoot  and  undershoot  of  approxi¬ 
mately  60  percent,  the  frequency  of  oscillation  is  slower  than  that  of  ,5 
second  time  delayed  system;  Fig,  8.7(d)  shows  the  plant  output,  rate  out¬ 
put  and  system  error  for  a  system  with  delay  time  of  2,5  seconds  and  k^_  = 
.02,  the  plant  output  continuously  oscillates  with  an  overshoot  and  under¬ 
shoot  of  approximately  135  percent,  the  frequency  of  oscillations  is  much 
slower  than  that  of  .5  second  time  delayed  system, 

8.5  Time  delay  in  the  phase  plane  diagram. 

The  effect  of  time  delay  in  application  or  removal  of  the  applied 
torque  after  the  control  has  signaled  that  the  operation  should  be  perform¬ 
ed.  It  will  be  of  interest  to  study  the  effect  in  terms  of  the  phase  plane 
diagram. 
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Fife.  8-6_ _ J  Time  rise  vs.  delay  time 


i 
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plant  output, velocity /Snd  system  error  of 'an  undeleyed  relay  servo  system 
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Fig.  8-7 b  •  The  plant  output, velocity  end  error  of  relay  servo  system  vdth  a  time  delay  of 
second  and  k  -  .02. 
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'c.  The  plant  output, velocity  and  system  error  of  a  relay  servo  system  with  a  tifl9  delay 
second  and  k  -  .02 


7 
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Because  of  the  time  delay  the  trajectory  in  the  phase  space  will  be 
prolonged  beyond  the  point  c°  =  R*  by  an  amount  such  that  the  difference  in 
coordinates  is  just  equal  to  the  time  delay*  that  in  Fig*  8*8  the  torque 

should  reverse  at  as  but  it  does  not  reverse  until  B  is  reached*  and  (C^  *  C^) 
(v®  -  v')/2  is  equal  to  t®  *  Similarly*,  the  torque  is  reversed  at  point  D 

D  n  Q 

instead  of  at  c5  and  at  F  instead  of  at  E*  It  is  apparent  that  a  curve  could 
be  drawn  through  points  such  as  B  and  F  so  that  the  switching  of  torque 
direction  would  take  place  when  the  trajectory  crossed  this  line8  and  the 
line  from  B  to  D  would  play  the  same  role  as  the  line  through  R°  for  the 
control  without  time  delay* 


Fig*  8*8.  Normalized  phase  plane  diagram  for  relay  servo 
system  with  a  time  delay  of  *1  seconds  for  a 
step  unit* 

The  effect  of  time  delay  is  thus  to  skew  the  switching  line  from  the 

verticle  in  the  phase  space s  and  quite  approximately  the  angle  of  skew  is 

-tan  and  is  in  clockwise  direction* 
a 

It  is  seen  in  Fig.  8*8  that  the  trajectory  appears  to  repeat  itself^ 
since  the  last  two  oscillations  are  nearly  identical*  The  control  oscillates 
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with  a  constant  amplitude  on  the  trajectory.  Such  a  repet ive  diagram  re¬ 
presenting  the  steady  state  oscillation  is  called  limit  cycle.  From  the 


experimental  work  it  shows  that  the  relay  controls  with  time  delay  and  nr. 
dead  zone  exhibits  the  defect  of  continuous  oscillate  loins . 

In  Fig*  8.9  is  shown  a  family  of  trajectories  so  arranged  that  the 
initial  and  final  points  of  each  are  symmetrically  located  with  respect  to 
the  origin.  Thus  each  trajectory  represents  just  such  a  limit  cycle*  and 
since  the  time  delay  for  each  is  fixed*  the  locus  of  the  end  points  of  these 
limit  cycles  measures  the  time  delay  required  to  obtain  the  amplitude  of 
oscillation.  As  would  be  expected*  the  larger  the  time  delay  the  larger  is 
the  magnitude  of  the  oscillation*  and  with  that*  the  lower  is  the  frequency 
of  oscillation  as  indicated  in  Fig,  8„5„  In  Fig.  8.10  is  shown  a  plot  of 
time  delay  and  oscillation  amplitude 0  It  turns  out  to  be  a  straight  line* 
therefore,  ve  can  obtain  an  equation  to  express  the  relation  of  oscillation 
amplitude  with  time  delay  as 


A  esc  -  *  & 


(8  A) 


the  above  equation  is  for  a  second  order  system  specified  by 


It  Is  interesting  to  see  the  relay  type  control  that  can  continuosuly 
cycle  when  a  time  delay  is  present.  The  relay  type  of  control  always 
oscillates  about  the  equilibrium  position.  For  any  time  delay  the  system 
would  always  oscillate*  and  with  a  magnitude  of  oscillation  depending  on 
the  time  delay.  The  relay  type  control  can  never  be  unstable  in  the  sense 
that  the  amplitude  of  oscillation  will  continue  to  increase  without  limit. 
It  is  obvious  that  this  is  the  advantage  of  the  relay  servo  system. 

In  Fig,  8.11(a)  the  phase  plane  of  the  relay  servo  system  with  a  time 
delay  of  zero  and  .5  seconds  is  plotted  on  the  same  sheet.  In  Fig,  8, 11(b) 
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Fig.'S -9.  Phase  plane  trajectory  showing  loci  of  limit  cycle  for  oscillation  of  constant 
amplitude  with  various  time  delay 
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Fig.  8  10c  Curve  showing  the  relation  between  time  delay  8nd  amplitude  of  oscillation 
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Fig.'  8-lla  •  Phase  plane  diagram  for  e  relay  servo  ssystem  showing  the  effect  of  time  delay 


Fig.  E-lib.  Phase  plane  diagram  for  a  relay  servo  system  showing  the  effect  of  time 
delay 


is  the  phase  plane  of  the  relay  servo  system  a  time  delay  of  one  and  2*5 
seconds  plotted  on  the  same  sheet.  From  these  phase  portraits  we  find  a 
most  important  aspect  of  the  phase  portrait  is  the  switching  line  which  is 
not  vertical.  The  effect  of  the  dead  time  is  to  rotate  the  switching  line 
into  the  first  and  third  quadrants*  the  result  is  growing  oscillation* 

The  process  of  moving  the  switching  line  from  vertical  axis  to  a  line 
in  the  first  and  third  quadrants  is  seen  to  affect  the  nature  of  the  response 
very  unfavorably.  On  the  other  hand*  if  the  switching  line  can  be  made  to 
lie  in  the  second  and  fourth  quadrants  by  means  of  phase  advance*  or  lead., 
on  the  error  signal*  the  relay  servo  system  is  likely  to  be  a  very  well  be¬ 
haved  one*  in  that  it  has  small  overshoot  and  will  reach  and  remain  within 
a  given  value  of  error  in  a  short  time.  This  will  be  seen  in  section  9* 

Dead  time  will  result  in  a  terminal  limit  cycle  which  may  or  may  not  be 
acceptably  small*  the  amplitude  of  oscillation  depends  upon  delay  time  as 
shown  in  Fig.  8.9  and  .10* 

8.6.  The  effect  of  time  delay  in  steady  state  oscillation  frequency. 

In  Fig.  8.12  the  curve  of  system  error  against  time  is  shown*  Initial 
condition  of  the  error  is  due  to  a  suddenly  applied  step  input  signal*  this 
gives  an  initial  error  and  velocity  (output  rate)  to  the  form  of  response 
waveform  as  shown  in  Fig.  8.12.  After  the  transient  has  died  out  the  out¬ 
put  shaft  follows  the  input  by  performing  oscillation  about  a  final  posi¬ 
tion  as  shown  in  Fig.  8.5. 

In  Fig.  8.12  curves  1*  2*  3,  4*  and  5  are  for  a  delay  time  of  .1*  -25* 

.5,  .8  and  1.0  second  respectively.  It  indicates  that  the  smaller  time 
delay  causes  high  frequency  oscillation  with  smaller  amplitude*  in  other 
words*  the  frequency  of  the  oscillation  is  higher*  and  the  larger  value 
of  time  delay  results  in  larger  oscillation  amplitude  with  lower  frequency 
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as  shown  in  Fig,  8.12, 

In  Fig.  8.13  curves  1,  2,  3,  4,  5,  and  6  are  for  a  delay  time  of  0>  .  1, 
.25*  ,5,  .8,  and  1.0  second  respectively.  The  curves  show  the  relation  of 
the  rate  output  to  time  delay. 

From  Fig.  8.5  or  Fig.  8,12S  the  data  for  oscillation  frequency  against 
time  delay  can  be  obtained  as  shown  in  Table  8.2. 

Table  8.2,  Fata  for  steady  state  oscillation  frequency  vs.  time  delay. 


delay  time 

.1 

.25 

.5 

.8 

1.0 


frequency 

,455 

,294 

.200 

.156 

.135 


The  data  plotted  in  log  log  scale  sheet  there  appears  to  be  a  straight 
line  as  shown  in  Fig.  8,14;  it  is  seen  that  the  frequency  decreases  as 
the  delay  time  increases.  Because  the  curve  of  oscillation  frequency 
against  time  delay  is  a  straight  line,  therefore  we  can  obtain  an  equation 
to  express  their  relation  as 


-/as 

f~-s7  (8.5) 

Ld 

From  Fig.  8.12  we  take  the  steady  state  error  for  each  time  delay, 
tabulate  the  data  as  shown  in  Table  8.3. 


Table  8.3.  Steady  state  error  vs.  time  delay. 


Delay  time  Error 

0  0 
.1  .09 

.  25  . 18 

.5  .31 

.8  . 484 

1.0  .588 


Fig.  8.15  is  the  plot  of  error  vs.  delay  time.  It  comes  out  approximately 
a  straight  line,  it  indicates  that  the  steady  state  error  is  proportional 
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l2JSrror  vs.  time  curves  showing  the  effect  of  time  delay  on  the  system  error 


^6 


Figj  8™ 13c Velocity  vs. time  curves  showing  the  effect  of  time  delay  on  velocity. 


Fig  8-14*.  Steady  state  oscillation 
frequency  vs.  time  delay. 


Fig. '8-15?  Systemerror  against  time  delay 


1.0 


Co  the  t i iK«ii  delay. 

8.7.  The  effect  of  the  magnitude  of  the  relay  voltage. 

From  the  experimental  work,  we  find  chat  the  magnitude  of  the  relay 
voltage  also  has  some  effects  on  the  response  curve  and  phase  trajectory  of 
the  system.  VJhen  the  relay  voltage  is  smaller  the  response  has  a  longer 
rise  time  and  a  smaller  overshoot  and  undershoot,  and  there  is  presented 
some  velocity  saturation  in  the  phase  trajectory*  the  amplitude  of  the 
limit  cycle  is  smaller. 

In  Fig.  8,16  shows  two  curves  with  the  same  delay  time  of  *5  second 
and  same  value  of  k,  -  ,02*  but  with  different  relay  voltage,  curve  (1)  is 
with  V  volt  and  curve  (2)  is  with  V  ~«5  volt.  From  these  curves  it  is 
indicated  that  the  smaller  relay  voltage  causes  smaller  oscillation  amplitude 
and  lower  frequency,  but  it  presents  a  longer  time  rise  which  is  its  dis¬ 
advantage;  so  far  as  the  time  rise  is  concerned  the  larger  value  of  relay 
voltage  is  desirable,  therefore,  in  all  the  experiments  we  set  the  relay 
voltage  equal  to  .5  volts.  The  rise  time  for  .2  volts  of  relay  output  is 
6.5  seconds,  and  the  rise  time  for  .5  volts  is  3.35  seconds.  The  overshoot 
and  undershoot  for  .2  volt  system  is  12,5  percent,  and  that  for  .5  volt  is 
30  percent*. 

In  Fig.  8.17  the  phase  plane  with  the  same  delay  time  and  same  value  of 
is  shown,  but  with  different  relay  voltage,  it  presents  a  velocity  satura¬ 
tion  for  the  system  with  .2  volts  of  relay  output,  and  the  limit  cycle  is 
smaller. 
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16. Curves  showing  the  effect  o£  the  relay  voltage  to  the  response.  t^=.5,k^=.02 
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Fig.  8-I?.Fhase  plane  diagram  showing  the  effect  of  relay  voltage,  t.r.5  second, kt=. 02. 


9 .  Elimination  of  the  steady  state  oscillation 


In  many  applications  steady  state  oscillations  of  the  output  shaft  are 
inadmissible;  besides  causing  continual  operation  of  relay  contacts  and 
wear  of  gears  and  bearings  the  power  consumed  in  oscillations  may  be  objec¬ 
tionable.  The  simplest  method  of  eliminating  the  oscillation  is  to  intro¬ 
duce  a  dead  zone  in  the  relays  so  that  the  error  always  returns  to  the 
dead  zone  after  a  disturbance.  But  the  relay  used  has  been  thus  far  assumed 
an  ideal  relay  with  no  dead  zone,  so  we  shall  not  discuss  this  possibility, 

9.1  Velocity  feedback  stabilization 

A  very  effective  method  of  reducing  the  settling  time  and  limit  cycle 
oscillation  of  an  relay-servo  is  to  make  the  relay  switching  conditional 
not  only  on  misalighment  but  also  on  its  rate  change.  This  can  be  achieved 
by  feedback  within  the  position  control  loop  if  the  feedback  signals  can 
be  added  algebraically  to  the  misalignment  signal  at  the  error  sensing 
device . 

The  effect  of  negative  velocity  feedback  is  illustrated  on  the  phase 
plane  diagram  in  Fig.  9.1;  and  the  corresponding  time  response  is  shown 
in  Fig.  9.2. 

It  is  obvious  that  the  larger  the  value  of  k  the  smaller  the  ampli¬ 
tude  of  the  limit  cycle,  and  the  more  the  switching  line  shifted  to  counter¬ 
clockwise  direction;  and  the  amplitude  of  oscillation  is  reduced;  the  over¬ 
shoot  and  undershoot  for  k^_  =  .02  is  about  15  percent,  for  k^  =  .8  is 
about  2  percent  as  shown  in  Fig.  9.2, 

9.2  Optimising  the  transient  response. 

Because  of  the  introduction  of  the  time  delay  in  the  system,  it  is 
difficult  to  derive  a  method  to  determine  the  value  of  k^  which  will  give 
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Fig.  9  =  1 0  Phase  plane  of  relay  servo  system  with  a  time  delay, of  .25  second  compensated 
by  increasing  the  gain  of  tachometer  feedback  gain  to  .8 
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Fig.  3-2.  Curve  showing  the  response  of  the  relay  servo  system  be  compensated  to  an  acceptable 
optimum  operation  by  in  creasing  kt  to  a  value  of  .8,  time  delay  .25  second. 


an  optimum  operation,  which  can  only  be  determined  through  experimental ion* 
As  we  know  the  system  damping  will  be  improved  by  the  tachometer  feed¬ 
back,  the  net  effect  of  the  derivative  signal  is  to  cause  the  relay  to 
operate  early,  so  by  using  k^  we  can  control  the  switching  characteristics 
of  the  relay.  For  many  applications  of  relay  servomechanism,  the  tacho¬ 
meter  feedback  compensation  is  adequate* 

Through  experimentation  it  was  found  that  the  system  with  a  delay  time 
below  ,5  second  can  be  compensated  for  5  percent  overshoot  and  undershoot 
which  is  considered  to  be  acceptable. 

The  response  and  phase  plane  of  the  compensated  system  for  various  time 
delays  are  shown  in  Fig*  9.3(a)  through  Fig*  9. 3(e)*  They  have  been  compen¬ 
sated  by  velocity  feedback.  The  values  of  k  (which  give  acceptable  optimum 

t 

operation)  are  as  follows: 

Table  9*1  Optimum  value  of  k  for  eliminating  the  effect 

of  time  delay* 

t  k  o/o  overshoot 

and  undershoot 


0 

.3 

/ 

.1 

.4 

.1 

.25 

.8 

1.6 

.5 

1.5 

5.2 

For  systems  with  .8  and  1.0  second  time  delay ,  it  is  found  that  they  can 

not  be  compensated  for  5  percent  overshoot  and  undershoot,  even  when  we 

increase  the  tachometer  gain,  ,  to  a  considerable  value,  the  amplitude 

of  steady  state  oscillation  is  still  rather  large  as  shown  in  Fig.  9*4  and 

Fig.  9*5;  the  only  effect  of  increasing  k  at  this  time  is  to  shift  the 

t 

whole  response  curve  downward,  as  shown  in  Fig*  9*5  the  response  of  the 

system  with  one  second  time  delay  at  k  -  4*5  is  rather  lew  but  still  with 

t 
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considerable  oscillation  amplitude  which  is  almost  the  same  as  that  when 
=  1.5,  which  is  beyond  the  acceptable  value. 

Fig.  9.6  illustrates  the  time  delay  vs,  curve  at  the  range  from  0  to 
.5  second  time  delay.  The  curve  appears  to  be  a  straight  line,  but  when 
the  time  delay  above  ,5  second  occurs  this  situation  does  not  hold,  since 
the  system  with  a  delay  time  above  .5  second  cannot  be  compensated.  When 
the  time  delay  is  less  than  .5  seconds,  the  system  can  be  compensated;  and 
because  t^  vs.  k^  curve  is  a  straight  line,  therefore,  it  can  be  expressed 
by  an  equation  as  given  in  equation  9.1. 

Kt=  -+  ^-8  td  (9.1) 

Table  9.2  and  9.3  are  the  results  for  comparing  the  uncompensated  and 
compensated  system.  The  data  in  Table  9.2  is  for  the  undelayed  system 
for  which  k^  =  . 3  will  compensate  the  system  to  give  optimum  operation. 

The  data  in  Table  9.3  is  for  a  .5  second  time  delayed  system  for  which  k^ 
=1.5  will  compensate  the  system  to  give  acceptable  optimum  operation. 
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plant  output , veloci ty  and  error  vs.  time 


Phase  plane 

Fig."9-3b."J_  Response  and  phase  plane,  t^  -  .1  sec. 
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plant  output  veloci ty  and  error  vs.  time 


phase  .plane 
Fig.  9- 3c.  4 


Response  and  phase  plane, t^  -  .25 
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plant  output, velocity  and  error  vs.  time 


Fig; J  REsponse  and  phase  plane, t^-  .5  sec.  -  1.5 
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Fig/  9-6.Tims  delay  vs.  tachbmoter  gain  curve 


I'afc  9,2  alculated  response  of  relay  servo  system  Forcing  function  = 
it  step  t  =  C  second 


T  itns 

System 

error 

Rate 

error 

Plant 

output 

k  =  .02 

t 

k  =  .  3 
t 

k  =  .02 

K  -  -3 

k  =  .02 

t 

k  -  . 
t 

0 

1 . 0000 

1,0000 

.0000 

.  0000 

.0000 

.  0000 

.  3 

.9467 

.9467 

-.1967 

- . 1967 

.0531 

.0531 

1.0 

.5161 

.8161 

-.3161 

-.3161 

.1839 

.1839 

1,3 

.5384 

.6384 

-.3884 

- . 3384 

.3616 

.3616 

2*0 

.4  323 

.4323 

-.4323 

-.4323 

.  5677 

.5677 

'N 

Z  0 

2090 

.  2090 

- . 4590 

-.4590 

.7910 

.7910 

3.0 

-.0223 

.0296 

-.4052 

-.1802 
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1C.  The  response  of  relay  servo  system  to  ramp  function  input. 


In  the  previous  sections  we  have  been  concerned  primarily  with  the 
response  to  step  function  input.  Experimental  results  obtained  from  the 
relay  servo  system  with  a  ramp  function  input  are  shown  in  several  graphs. 

When  a  unit  ramp  input  is  applied  to  the  relay  servo  system,  the  out¬ 
put  and  the  error  of  the  system  are  shown  as  a  straight  line,  and  the  velocity 
for  most  part  is  a  horizontal  line,  as  shown  in  Fig.  10.1.  The  Cdot  vs.  C 

and  Edot  vs.  E  phase  planes  have  the  same  shape,  merely  the  Cdot  vs.  C  plane 

appears  above  X-  axis.  Fig.  10.2  shows  the  Edot  vs.  E  plane.  This  is  with 
a  relay  voltage  equal  to  .5  volt.  When  the  relay  voltage  is  increased  to 
1.0  volt  the  same  situation  remains.  For  relay  voltage  at  .5  and  1.0  volt 

the  only  effect  of  the  time  delay  on  the  relay  servo  system  is  to  cause  some 

delay  in  the  response,  no  oscillation  has  been  observed. 

But  when  the  relay  voltage  is  increased  to  1.5  volts,  this  situation  is 

* 

changed.  For  an  undelayed  system  the  output  still  appears  as  a  straight  line 
but  the  system  error  is  no  more  a  straight  line,  and  the  velocity  curve  is 
with  some  oscillations  as  shown  in  Fig.  10.3.  Cdot  vs.  C  and  Edot  vs.  E 
phase  planes  have  a  different  shape  as  shown  in  Fig.  10.4  and  Fig.  10.5. 

For  delayed  system  the  output  no  more  remains  a  straight  line  but  with  some 
oscillations,  and  the  system  error  and  velocity  are  much  different  than  that 
of  the  system  with  relay  voltage  of  .5  or  1.0  volt  as  shown  in  Fig.  10.6. 

The  Cdot  vs.  C  and  Edot  vs.  E  phase  planes  are  quite  different  as  shown  in 
Fig.  10.7  and  Fig.  10.8  respectively.  Oscillations  and  limit  cycle  has  been 
observed  when  the  relay  voltage  is  1,5  volts. 

The  effect  of  time  delay  on  the  response  of  the  system  can  also  be  com¬ 
pensated  by  velocity  feedback.  For  the  .5  second  time  delayed  system,  when 
k  =  2.0,  the  response  can  be  compensated  to  an  acceptable  optimum  operation 
as  shown  in  Figures  10.9,10,  II. 
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Fig.  10-1.  C,  CD0T,E  vs.  time 
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Fig.  1 .  Phase  trajectory  of  relay  servo  system,  relay  voltage 
.5  volt  t- 
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Fig  15-7.  phase  plane  (rate  output  vs.  output)  of  relay  servo  system  with  unit  ramp  input, 
time  delay  .5  second, k  .02,  relay  voltage/. 5  volt. 


Fig.  !£■ Phase  trajectory  of  ralay  servo  system  to  an  unit 
ramp  input,  time  deiay  —  .5  second,  relay  voltage—/.  5  volt. 
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11.  Conclusions, 


The  aim  of  this  thesis  was  to  demonstrate  the  effects  of  transportation 
lag  on  the  relay  servo  system. 

For  conducting  the  experiment  by  means  of  an  analog  computer*  first  one 
could  simulate  the  time  delay  by  Pade*s  approximation  or  by  magnetic  tape 
recorder.  It  was  seen  that  the  generation  of  time  delay  from  magnetic  tape 
recorder  presents  a  more  accurate  and  better  result  than  that  from  the  Pad eBs 
approximation.  Therefore,  the  magnetic  tape  recorder  was  selected  to 
generate  the  time  delay  throughout  the  experiment. 

In  this  investigation  most  of  the  experiments  were  conducted  by  using 
a  digital  computer  which  gave  more  accurate  results  than  the  analog  com¬ 
puter  does;  especially  since  it  can  give  an  exact  value  at  any  desired  in¬ 
stant.  The  analog  computer  has  also  been  used  to  check  the  results  ob¬ 
tained  from  the  digital  computer.  It  is  noticed  that  they  are  quite  in 
agreement . 

From  the  analysis  and  experiments  the  following  conclusions  has  been 
obtained. 

The  root  locus  of  the  system  was  greatly  changed  due  to  the  intro¬ 
duction  of  the  transportation  lag;  and  because  time  delay  introduced 
periodicaly  a  phase  shift  of  360  degrees,  it  results  in  an  additional  (upper) 
branches  when  at  high  frequency. 

The  effect  of  time  delay  on  the  transient  response  is  to  increase  the 
overshoot  and  undershoot  and  causes  oscillations  and  also  leads  to  a  longer 
rise  time,  the  response  of  the  system  can  never  be  settled. 

The  effect  of  time  delay  on  the  phase  plane  is  to  rotate  the  switching 
line  into  the  first  and  third  quadrants,  the  results  is  growing  oscillations 
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In  other  words  the  transportation  lag  will  result  in  a  terminal  limit 
cycle  which  may  or  may  not  be  acceptably  small. 

The  plot  of  time  delay  versus  oscillation  amplitude  came  out  to  be  a 
straight  line;  the  delay  time  multiplied  by  a  constant  of  .6  gives  the 
magnitude  of  the  oscillation  amplitude. 

The  time  delay  versus  steady  state  oscillation  frequency  when  plotted 
on  log  log  coordinates  came  out  to  be  a  straight  line  in  such  a  relation 
that  a  constant  of  ,135  divided  by  gives  the  oscillation  frequency, 

A 

The  steady  state  error  and  the  velocity  are  proportional  to  the  time 
delay. 

These  effects  caused  by  the  transportation  lag  can  be  compensated  to 
some  extent.  The  method  used  in  this  investigation  to  eliminate  steady 
state  oscillation  and  to  optimise  the  transient  response  by  employing  tacho¬ 
meter  feedback.  It  was  seen  that  a  system  with  a  delay  time  below  ,5  second 
can  be  compensated  to  an  acceptable  small  oscillation  amplitude  (5  percent). 
But  for  time  delay  above  .8  second  it  cannot  be  compensated  to  such  a  small 
amount;  the  only  effect  of  increasing  tachometer  gain  is  to  shift  the  whole 
response  curve  downward  with  the  same  amount  of  steady  state  over-  and  under¬ 
shoot  . 

The  system  with  a  time  delay  below  .5  second  can  be  compensated  to  reach 
an  acceptable  optimum  operation.  For  delay  time  below  0 5  second  its  relation 
with  tachometer  gain  is  approximately  a  straight  line,  therefore,  it  can  be 
expressed  by  an  equation. 

When  a  time  delay  is  present  the  relay  type  control  can  continuously 
oscillate,  and  with  a  magnitude  of  oscillation  depending  on  the  time  delay. 
The  relay  type  control  can  never  be  unstable  in  the  sense  that  the  amplitude 
of  oscillation  will  continue  to  increase  without  limit.  This  is  one  advan¬ 
tage  of  the  relay  servo  system. 
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